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ABSTRACT 

The  measurement  of  residual  stress  by  the  X-ray  diffraction  technique 
involves  determining  the  lattice  spacing  of  a crystallographic  plane  at 
different  Inclinations  of  the  sample  and  relating  the  change  in  the  spacing 
to  a stress  on  the  surface  of  the  specimen.  The  sources  of  error  in 
determining  the  residual  stress  are  investigated  in  this  report.  A short 
review  of  the  fundamental  principles  of  the  measurement  is  given.  The 
Important  instrumental  and  geometric  factors  contributing  to  errors  in 
the  measured  stress  are  presented.  To  account  for  random  errors  in  the 
data  accumulation  on  the  measured  stress,  a complete  statistical  analysis 
based  on  a least-squares  parabola  is  given. 

Employing  an  automated  diffractometer  (see  T.R.  No.  16),  an  extensive 

investigation  of  the  precision  of  the  residual  stress  measurement  was 

made  to  evaluate  the  various  procedures  offered  in  the  literature  to 

optimize  the  automation.  The  results  Indicate  that  a multiple  least- 

squares  parabola  is  the  most  reproducible  method  of  defining  the  peak  of 

the  profile,  in  contrast  to  the  common  three  point  parabola  used  in  the 

United  States,  the  center  of  gravity  used  in  Germany  or  the  middle  of  the 

half-maximum  intensity  used  in  Japan.  Also,  for  a standard  diffractometer, 
2 

the  sin  f technique  is  more  precise  than  the  usual  two  tilt  method  even 


I 
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when  the  total  time  of  analyals  for  each  is  identical.  In  studying 
various  beam  optic  arrangements,  it  was  found  that  the  stationary  slit 
(non-focusing  method)  offers  the  best  system  in  that  the  measured  stress 
is  only  mildly  sensitive  to  sample  displacement  (~  7 MPa  for  .25  non  dis- 
placement) and  introduces  very  little  Instrumental  broadening.  The 
operational  principles  and  use  of  a position  sensitive  detector  (PSD) 
for  the  measurement  of  residual  stress  has  been  described  previously  (see 
T.R.  No.  11).  The  PSD  simultaneously  collects  data  over  the  entire 
diffraction  profile  enabling  the  use  of  a multiple  point  parabola  to 
define  the  peak  without  having  to  step  scan  the  profile.  Extensive 
testing  is  described  which  showed  that  this  detector  could  achieve  a 
speed  of  analysis  heretofore  never  obtainable  in  the  measurement  of 
residual  stress. 
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CHAPTER  1 

The  X-ray  diffraction  procedure  for  determining  the  surface  residual 
stresses  Is  well  established.^-7^  Each  family  of  Identical  planes  of 
atoms  In  a polycrystalline  material  has  a constant  Interplanar  spacing. 


dhki  whlch»  *h*n  act«d  upon  by  an  elastic  stress,  changes  to  a new  value 
dependent  on  the  direction  and  magnitude  of  that  stress.  A change, 

Adhki,  in  the  interplanar  spacing  will  cause  a corresponding  change, 

&0,  in  the  Bragg  angle  of  diffraction  by  the  family  of  planes.  The  strain, 
Ad/d,  can  be  measured  by  the  change  in  the  diffraction  angle  and  the 
stress  can  be  obtained  from  the  strain  with  formulae  derived  from  linear 
Isotropic  elasticity  theory. 

The  principal  stresses  9.  and  o0,  the  surface  stress,  o , and  the 

1 « cp 

corresponding  strain  are  shown  In  Fig.  1.1.  The  term  f is  the  angle 
between  the  surface  normal  H the  direction  of  strain  being  measured. 
Application  of  isotropic  elasticity  theory  yields  i>2  following  relation- 
ship between  the  principle  stresses,  surface  stress  and  measured  strain. 


i*vl°2*  -1  fa, 


(1.3) 


▼s  this  equation,  d Is  the  lattice  spacing  In  the  direction  defined  by 
<P»f 

9 and  t (see  Fig.  1.1)  and  dQ  la  the  interplanar  spacing  of  the  stress 

free  state. 

Equation  1.3  forma  the  basis  of  X-ray  stress  analysis  and  is  utilised 
primarily  in  two  methods;  the  'sin2^'  and  the  'two  tilt'  technique. 
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(A)  'Sin2^'  Technique 

In  this  method  several  value*  of  lattice  strain  are  measured,  each 

at  a different  f tilt  of  the  specimen.  It  is  possible  to  determine  any 

surface  component  of  stress  from  a least-squares  straight  line  for  the 

2 

lattice  strain  as  a function  of  sin  if. 

From  Eq.  1.3: 


3 «in  if 


(1.4) 


In  terms  of  the  interplanar  spacing  it  follows  that: 


3 sin2f 


(1.5) 


Because  several  values  of  d are  determined,  errors  resulting  from 

9 

random  fluctuations  of  d are  minimised.  Four  or  si*  * tilts,  taken 

*P»? 

in  equal  increments  of  sin  f,  are  normally  utilized. 

(B)  'TWo  Tilt*  Technique 

Isotropic  elasticity  theory  predicts  the  strain  < to  be  linearly 

2 ®’f 

dependent  on  sin  f,  as  in  Eq.  1.3.  When  this  holds  true  only  two  incli- 
nations of  the  sample  are  necessery  to  determine  the  surface  stress.  The 
Interplanar  spacing*  are  determined  at  f-0  and  at  an  inclination  of 
The  formula  relating  the  stress  to  the  strain  is  than  given  by: ^ 


d_  -d 


(i.«) 
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The  c*ria  (E/l+v) • l/«in  f is  often  combined  Into  a calibration  constant, 
K,  which  can  be  experimentally  determined  for  a particular  combination 
of  f and  reflecting  planes  in  a given  material.  Experimental  determination 
of  K Is  desirable  because  bulk  values  of  E and  v are  not  necessarily 
applicable. ^ 

It  Is  unfortunately  cooeon  practice  to  replace  (d  -d  )/d  in  Eq.  1.6 

tp.f  o o 

by  the  approximation  -cot0 * ^(29^20^)  on  the  basis  of  Bragg's  law  to 
obtain  a formula  in  terms  of  the  peak  position  20: 


% 


tt  r e 

no  • kjz)  * 


cot  *<9o+0t) • (29o-29+> 


(1.7) 


where  20  and  29  ere  In  degrees.  For  ssmII  peak  shifts  this  does  not 

o * 

Induce  much  error.  The  stress  constant  then  becomes: 


„ 1 . TT 

K"  2 180  * 


cot  %(0o+9j) 


(1.8) 


and 


a - K A20  . 

<P 

CHAPTER  2 

The  methods  of  residual  stress  measurement  by  X-ray  diffraction 
just  described  are  based  on  the  fundamental  equation  (Eq.  1.3)  derived 
by  Isotropic  elasticity  theory  and  assume  homogeneous  deformation.  These 
methods  have  been  shown  to  sometimes  yield  anomalous  results  when  applied 
to  samples  which  have  been  plastically  deformed  unlaxlally  in  tension^18. 
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n.9  20} 

or  by  rolling.  ’ ' The  anomalous  results  ere  not  true  mechanical 

macrostresses  and  may  be  attributed  to  a number  of  causes. 

nq\ 

Recent  work  by  Marion  and  Cohen v ' has  led  to  the  quantification 
of  a deformation  model  to  account  for  one  of  the  anomalies  in  the  X-ray 
technique  of  residual  stress  measurement  which  occurs  after  plastic 
deformation.  The  classical  formulae  around  which  all  X-ray  technique  of 

residual  stress  analysis  (Eq.  1.3)  predicts  a linear  dependence  of  d 



2 

on  *ln  f . A non- linear  dependence  of  the  X-r*y  measured  strain  on 
2 

sin  f has  been  reported  for  a number  of  materials  which  have  undergone 
plastic  deformation  due  to  elongation^*  ^ or  rolling.  ^*®^  Such 
deviations  from  the  theory  has  prompted  questions  concerning  the 

(19) 

validity  of  the  X-ray  measurement  in  such  cases.  Marion  and  Cohen 

2 

attributed  the  non-linear  dependence  of  d on  sin  f to  the  relief  of 

microstrelns  in  subgrain  Interiors  which  are  oriented  to  be  relieved  by 

(23  24) 

a dynamic  recovery  process  proposed  by  Veldemann  et  al.  ’ ' This 

produces  a non-random  distribution  in  the  interplanar  spacing  which  is 

related  to  the  texture  developed  during  the  plastic  deformation  process. 

The  authors  developed  a distribution  function,  f(f),  describing  the 

variation  in  interplanar  spacing  at  each  f inclination,  and  incorporated 

2 

the  non-linear  behavior  of  d vs.  sin  f into  the  general  formulae  for  the 

X-ray  method  of  residual  stress  analysis.  By  measuring  both  the  inter- 

planar  spacing,  d , and  the  distribution  function,  f(f),  as  a function 
2 

of  sin  f,  the  non-linear  dependence  of  d may  be  separated  from  the  linear 
component  through  the  following  formula: 


v,  • <ws  >«♦>  ♦ 4»  Q&>  v1"**  * "•  • 


(2.1) 
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Th«  term  corresponds  to  the  lattice  spacing  in  a region  that  is 

fully  relieved  and  dj  the  lattice  spacing  in  a region  that  has  not  relieved. 

The  distribution  function  describes  the  variation  of  relief  with  orientation 

and  is  calculated  by  determining  the  texture  in  the  region  of  the  pole 

figure  for  which  the  residual  lattice  strain  Is  measured.  This  is  most 

easily  accomplished  by  measuring  the  Integrated  Intensity  of  the  diffraction 

peak  of  Interest  at  each  f inclination  and  normalizing  the  distribution 

function  by  setting  f(f)  » 1 at  the  maximum  value  of  the  curve  of 

integrated  intenalty.  The  correlation  between  the  change  in  the  distri- 

2 

butlon  function  and  the  osclllatlona  in  d vs.  sin  f la  shown  vividly  in 
Fig.  1.3. 
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FIGURE  1.3  * from  Morion'"' 

a)  d v#.  tin*  9 for  Armco  Iron  apaciman; 
raducad  69  pet  by  rolling;  211  peak 

. with  CtKq,. 

b)  Taxtura  for  sampla  daacribed  in  (a). 
f(t>  rapraaants  tha  norma li sad  taxtura 
distribution  function. 


1.1703 


1.1702 


1.1701 


1.1700 


1.1699 


1.1698 


10 


CHAPTER  3 

AN  EXAMINATION  OF  THE  ACCURACY  OF  X-RAY  RESIDUAL  STRESS  ANALYSIS 
3.1  INTRODUCTION 

The  confidence  In  any  measurement  depends  on  a critical  consideration 
of  the  sources  of  error  contributing  to  the  overall  accuracy.  The 
accuracy  may  be  separated  Into  two  components; 

1)  the  precision  or  closeness  of  agreement  among  repeated 
measurements,  and 

2)  the  bias  or  difference  between  the  mean  of  the  repeated 
measurements  and  the  true  value. 

Many  factors  are  known  to  contribute  to  the  bias  of  the  residual 
stress  measurement  by  X-ray  diffraction  techniques.  The  factors  are 
broadly  classified  by  Jatezak  and  Boehm  into  three  categories, 
equipment  and  Instrumental  factors,  technique  or  geometrical  factors  and 
specimen  factors.  The  most  important  Individual  errors  are  well  treated 
in  the  literature  and  a summary  is  presented  to  outline  the  source  and 
magnitude  of  each  contributing  factor.  Experimental  tests  are  described 
to  check  the  theoretical  treatments  of  the  most  important  error,  that  of 
sample  displacement.  These  treatments  account  for  correctable  errors  in 
the  measurement  but  do  not  delineate  the  potantial  accuracy  of  the  technique 
or  apparatus. 

The  precision  of  the  X-ray  residual  stress  measurement  relies  on 

reproducible  determination  of  the  peak  shift.  Numerical  methods  of  peak 

• OOia  3Wl  ttOi3vtt  H 1 r ; . 

location  are  mandatory  for  fully  coaputerlsed  control  of  the  measurement. 

The  peak  may  be  found  by  fitting  a smooth  curve  to  the  experimental  data. 
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Such  a curve  cannot  give  an  absolute  measure  of  the  peak  location  as  the 
data  is  subject  to  random  errors.  These  random  errors  occur  because  the 
diffracted  Intensity  is  measured  at  a finite  number  of  20  positions  and 
every  measurement  is  subject  to  a statistical  counting  error.  Statistical 
treatments  for  different  numerical  methods  for  measuring  line- prof lies 
can  be  found  in  the  literature.  These  include  the  centroid, 

(ftO 

median1'  * and  geometrical  peak.  The  peak  position,  located  by 

a least  squares  parabola  is  shown  here  to  be  the  most  precise  method.  A 
statistical  analysis  based  on  the  least-squares  parabola  is  extended  here 
to  the  determination  of  residual  stress.  The  region  of  fit  of  a parabola 
to  the  diffraction  profile  is  discussed.  The  commonly  used  three  point 
parabolic  fit  is  examined  theoretically  and  experimentally.  A test  on  the 
effect  of  the  number  of  data  points  is  made  to  determine  if  the  three 
point  fit  is  actually  the  optimum  procedure  for  recording  data.  It  is 
shown  that  for  fixed  total  time  of  data  accumulation,  the  observed  pre- 
cision from  replicate  tests  is  improved  by  using  many  data  points, 
especially  for  broad  profiles.  The  three  point  parabolic  fit  is  adequate 
only  for  samples  exhibiting  a sharp  profile. 

2 

The  precision  of  the  'two  tilt'  and  'sin  f'  stress  techniques  are 

tested  under  identical  conditions  so  that  comparisons  on  the  relative 

2 

accuracy  of  each  method  can  be  made.  The  'sin  f'  method  is  shown  to 
have  a greater  precision  that  the  'two  tilt'  procedure  when  a fixed 
total  time  of  data  accumulation  is  used.  In  addition,  the  precision  of 

★ 

Any  consistent  feature  of  the  line  profile  may  be  utilised  to  define 
the  reference  point  from  which  the  position  of  the  profile  is  based. 

Three  comaw n features  in  residual  stress  analysis  are  the  centroid, 
midpoint  of  the  half-intensity  breadth  value  and  the  peak.  The  posi- 
tions located  by  each  method  are  not  identical  but  must  be  reproducible. 
Comparisons  of  the  reproducibility  of  each  method  are  given  in  Sec.  3.3.2 
and  show  that  a least  squares  parabola  exhibits  the  best  reproducibility. 
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the  residual  stress  measurement  using  a position  sensitive  detector  is 
tested.  Because  many  data  points  are  obtained  in  the  region  of  curve 
fitting,  the  PSD  is  shown  to  be  a both  rapid  and  precise  means  of 
accumulating  the  data.  A least-squares  parabola  fit  to  the  region  of  the 
peak  of  the  diffraction  profile  obtained  with  the  PSD  is  shown  to  be  the 
most  reproducible  means  of  locating  a reference  position  on  the  profile. 

3.2  FACTORS  INFLUENCING  BIAS  IN  THE  MEASUREMENT 

There  are  three  broad  categories  into  which  all  factors  influencing 

the  bias  can  be  classified:  Instrumental,  geometrical,  and  specimen 

(29) 

factors.  A detailed  list  compiled  by  Jatezak  and  Boehm  is  reproduced 
in  Table  3.1. 

Many  of  the  sources  of  error  are  dependent  on  either  the  capability 
of  the  experimenter  (alignment,  sample  position)  or  operating  conditions 
(electronic  stability)  of  the  equipment  and,  as  suctv  are  unpredictable. 

An  explanation  of  the  predictable  factors,  their  effects  and  correction 
terms  are  given  below. 

3.2.1  Angular  Dependent  Intensity  Factors 

(A)  Lorentz- Polarization  Factor 

This  term  combines  two  20  dependent  intensity  corrections,  the 

lorentz  factor  and  the  polarization  factor,  which  arise  from  the  geometry 

(33) 

of  the  diffraction  process.  The  Lorentz  factor  for  a powder  can  be 

i 

conaldered  as  two  Independent  factors,  one  arising  from  the  number  of 
particles  in  the  sample  which  are  in  orientation  to  contribute  to  a 
particular  reflection  and  the  aecond  arising  from  the  fraction  of  the 
reflection  ring  which  is  detected.  The  first  factor  ia  independent 

At  , ilC  hftrffl  fKfM  I.  Vl  ' . ji  • T 'SI  1 T tiff.  * | 7 « 

f*  _ « • . ( •*  , •/  rrtif/j  ef,v  a fa-  , , , • , - r.  «u  * • «.  * , • 
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of  Che  detection  geometry  whereas  the  second  factor  is  dependent  on  the 

acceptance  function  of  the  detector and  the  velocity  of  the  reflection 

(33^ 

through  the  Ewald  sphere.  The  polarization  tens  arises  because  X-rays 

are  polarized  after  being  scattered  or  diffracted,  the  amount  of  polari- 
zation depending  upon  the  angle  through  which  it  is  scattered  or  diffracted. 
The  conventional  correction  factor  (combining  the  Lorentz  and  polarization 
factors)  for  filtered  radiation  and  point  counting,  designated  LP,  is:^7* 

LP  - -1-  ±-£.°.8/.2-0  . (3.1) 

sin  0 

The  measured  Intensity  need  only  be  divided  by  Eq.  3.1  prior  to  processing. 

Cooper  and  Glasspool^^  have  shown  that  in  some  cases,  notably 
when  the  axial  (perpendicular  to  the  plane  of  the  diffractometer)  half- 
height  of  the  detector  aperature  is  large  and  the  scattering  angle 
becomes  small  (0  -•  0°)  or  large  (0  - 90°),  the  curvature  of  the  dif- 

D O 

fraction  cone  becomes  Important.  In  the  traditional  derivation  of  the 
Lorentz  factor,  the  curvature  of  the  cone  is  Ignored.  For  20  values 
greater  than  about  9*  and  less  than  171°  the  conventional  Lorentz  factor 
is  in  error  by  less  than  1%,  as  calculated  by  Cooper  and  Glasspool. 

This  small  error  confined  with  the  fact  that  in  peak  location  measurements 
only  the  angular  dependent  changes  in  the  correction  factor  are  Important 
isake  the  conventional  Lorentz  factor  quite  satisfactory. 

(B)  Absorption  Factor 

The  path  length  of  the  primary  and  diffracted  X-rays  within  the 
specimen  differ  when  the  specimen  is  tilt  by  For  a flat  specimen  the 
measured  intensity  is  corrected  for  absorption  by  dividing  by:*  ’ 
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ABS  ■ 1 - tan  * cot  0 (3.2) 

Die  angle  f la  defined  aa  positive  when  the  normal  to  the  aample  moves 
towards  the  X-ray  aource.  Although  ^ la  a polar  angle,  lta  algn  muat 
be  known  In  order  to  calculate  the  correct  absorption  coefficient. 

(C)  Atomic  Scattering  Factor 

The  intensity  of  the  diffracted  beam  depends  in  part  on  the  square 

of  the  structure  factor.  The  angular  dependent  part  of  this  term,  denoted 

the  scattering  factor  f,  should  be  accounted  for  by  dividing  the  intensity 
2 

by  f . The  atomic  scattering  factor  is  a function  of  the  atoms  comprising 
the  crystal,  the  diffraction  angle  0 and  the  wavelength  of  the  radiation 
and  values  are  published  in  the  International  Tables  for  X-Ray  Crystallo- 
graphy. <65> 

(D)  Temperature  Factor 

The  Debye  temperature  factor  ia  the  reduction  in  intensity  dua  to 
thermal  motion  of  the  lattice.  Atoms  in  a crystal  vibrate  due  to  temperature 
so  that  at  any  Instant  corresponding  atoms  are  not  separated  by  exact 
multiples  of  cell  dimensions.  Scattering  from  two  such  atoms  will  not  be 
exactly  in  phase  resulting  in  a reduction  in  the  structure  factor.  The 
correction  for  this  reduction,  called  the  temperature  factor,  is  slightly 
dependent  on  sin  0 and  strictly  speaking  should  be  accounted  for. 

Over  the  range  of  154° 20  to  158° 20  where  the  meet  commonly  used 

analysis  of  ateel  occurs  with  Cr.  radiation.  Short  and 


~ ‘ — * 
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(28) 

Kelly'  calculated  the  change  In  the  intensity  factors.  Their  values 
are  reproduced  in  Table  3.2  and  indicate  that  the  Debye  temperature  factor 
is  very  small  in  relation  to  the  others.  For  generality,  the  Lorentz- 
po lari ration,  absorption  and  scattering  factors  have  been  Included  in 
the  automated  residual  stress  analysis  program  developed  in  this  study. 

3.2.2  Factors  for  "Beam  Optics" 

Beam  focusing  depends  on  the  geosmtric  arrangement  of  the  X-ray 
path  and  on  the  horizontal  and  vertical  divergence  of  the  beam.  The 
Important  smthods  of  'beam  optics'  are  discussed  below. 


(A)  Variation  of  the  Focal  Point  with  0 and  ^ 

There  are  two  methods  of  beam  focusing  used  in  X-ray  stress  measurement, 
the  para focusing  method  and  the  stationary  method.  An  excellent  study 
on  the  errors  associated  with  each  method  is  given  by  Zantopulos  and 


In  the  parafocus  technique  the  receiving  slit  and/or  detector  are 
moved  along  a radius  of  the  goniometer  toward  the  specimen  to  fulfill 
the  changing  beam  focusing  condition  when  the  inclination,  f,  of  the 
specimen  is  varied. *26) 


The  distance  from  the  sample  to  the  focus  position  is 

given  by: ^ 


r sssli.  +.  .Iff- 

QC  cos{f  - (90-0)) 


(3.3) 


where  is  the  goniometer  radius. 

In  the  stationary  method  the  receiving  slit  and  detector  remain  on 
the  goniometer  circle  at  all  times  thereby  deliberately  not  fulfilling 
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TABLE  3.2 

CHANGE  IN  CORRECTION  FACTORS  FROM  154°  TO  158° 20 
FOR  Fe  AND  CrK  RADIATION 


Factor 

154°  20 

158°  20 

Change 



2 

Lorentc  (1/sln  0) 

1.0533 

1.0378 

- 1.477. 

2 

Polarisation  (1+cos  20) 

1.8078 

1.8596 

+ 2.87% 

Absorption  (f  - (f) 

1. 

1. 

0 

Absorption  (f  - 45°) 

0.7691 

0.8056 

+ 4.74% 

Atoaic  Scat taring  (f2) 

1.3503 

1.3322 

- 1.34% 

Teaperature 

0.8653 

0.8634 

- 0.22% 

focusing  conditions.  A sacrifice  In  Intensity  Is  made  but  the  complication 
of  moving  the  receiving  slit  and/or  detector  la  avoided.  The  Japanese 
use  a version  of  this  technique  described  as  the  parallel  beam  method.  ' 
Using  long  Soller  baffles  or  plates  perpendicular  to  the  diffractometer 
plane  rather  than  in  the  usual  horltontal  position,  the  X-rays  are  made 
to  be  highly  collimated  parallel  beams  which  do  not  have  a focal  point. 

The  angle  of  diffraction  is  uniquely  defined  by  the  angle  between  the 
primary  and  diffracted  beam  (see  Fig.  3.1)  and  as  such  a receiving  slit 
is  not  necessary. 

(B)  Horizontal  Beam  Divergence 

In  the  parafocusing  method,  true  focus  demands  that  the  sample 
surface  lie  on  the  focusing  circle  which  is  given  by: ^ 

RpC  - Rqc/2  «io  (0  + *>  (3.4) 

True  focusing  demands  a continuous  change  in  the  curvature  of  the  apeclmen 
during  0 and  f angular  movements.  Since  this  is  generally  not  practical 
an  error  will  arise  which  is  dependent  on  the  curvature  of  the  sample 
and  the  horizontal  beam  divergence.  Marion has  derived  a simple 
formula  to  estimate  the  error  in  the  peak  shift  due  to  beam  divergence. 
Defining  a ••  half  the  angular  beam  divergence  the  peak  shift  in  degrees 
20,  6(A23)Bjj,  between  the  f ■ CP  and  t ■ t°  inclination  is: 
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FIGURE  3.1  Illustration  of  parallal  beam  geometry.  The 
angle  20  la  uniquely  defined  by  the  angular 
relationship  between  the  two  seta  of  aoller 
slits  and  is  independent  of  the  sample  position 
The  dashed  lines  represent  other  sample 
positions  that  would  not  affect  the  angle  of 
diffraction  . 


ao 


<>? 


— 

'*■  ! 


riur.  J20  . 02.  . . tMst  .cMtoy^-l  . 

t n R ' 2 Lco*(-o+f-cp)  cos  (cH-f-cp)  J Y 

The  Cera  R ' is  given  In  Eq.  3.3. 

P 

Because  there  is  e distribution  of  intensity  between  the  central 
beam  end  the  left  end  right  portions,  the  actual  peek  shift  will  be  less 
than  6(A20>Bq.  It  has  been  shown^*^  that  the  centroid  of  the  diffracted 
beam  approached  a limit  of  1/3  of  the  6(A20>bq  value  for  the  f,  a and  20 
values  used  in  Ref.  31.  (The  centroid  was  determined  by  a computer 
Iteration  method  in  which  the  divergence  angle,  cr*  was  divided  into  1000 
parts).  Therefore,  as  a conservative  estimate  one  can  use  1/2  of  the 
value  calculated  in  Eq.  3.5.  An  estimate  of  the  magnitude  of  this 
quantity  is  given  in  Table  3.2. 


(C)  Vertical  Beam  Divergence 

Grains  that  have  planes  slightly  tilted  from  that  of  the  diffracting 
position  for  a parallel  incident  beam  My  contribute  weakly  to  a peak 
giving  rise  to  an  apparent  peak  shift.  This  has  been  thoroughly  treated 
by  Cohen.  The  peak  shift  depends  on  the  amount  of  texturing  and  the 

slit  system  and  is  difficult  to  determine  exactly  but  is  quite  sasll. 

An  estimate  is  given  in  Table  3.2. 


In  the  parafocusing  end  stationary  slit  techniques  systeMtlc  errors 
are  produced  by  focusing  aberrations  during  diffraction  due  to  imperfect 
specimen  contours  and  beam  divergence  as  shown  in  Fig.  3.2.  The  photon 
rays  are  in  a divergent  beam  and  diffract  from  a sample  which  is  not 
curved  to  the  focusing  circle  so  that  the  point  of  focus  in  Fig.  3.2, 


TYPICAL  INSTRUMENTAL  ERRORS  IN  PEAK  LOCATION  FOR  20 
(Approximately  the  211  CrK  Peak  from  Iron) 
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Note:  Maximum  error  is  either  one  of  these  but  not  the  total  range. 

Calculated  for  steel  from  (20x  - 20  ) where  kA5  = 593  MPa/°2G  and  = 396  MPa/°20. 
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Is  not  really  a point  but  simply  a converging  area  for  the  beam  which 

(31) 

causes  a focusing  aberration.  Zantopulos  and  Jatczak  studied  these 
systematic  errors  by  analytically  describing  the  beam  path  for  any 
particular  ray  within  the  divergent  beam  from  the  source  to  the  detector. 

By  dividing  the  divergent  beam  into  1000  such  rays  they  calculated  the 
centroid  of  the  rays  at  the  detector  for  both  'beam  optic*  techniques 
using  specimens  of  varying  curvature.  They  concluded  the  focusing  error 
due  to  imperfect  specimen  contour  and  horizontal  beam  divergence  for  the 
stationary  slit  method  yields  about  one  third  the  error  of  the  para- 
focusing  method  at  20  - 156°;  the  discrepancy  increasing  for  smaller 
20.  in  addition,  the  error  when  using  a beam  divergence  of  1°  on  flat 
samples  or  samples  having  reasonable  curvature,  will  be  less  than  16  MPa 
(2.3  kail. 

One  factor  not  Included  in  the  analysis  by  Zantopulos  and  Jatczak 
is  the  ability  to  accurately  position  the  receiving  slit  when  the  para- 
focusing  technique  is  used.  The  exact  position  along  the  radius  is  not 
critical  because,  as  shown  above,  the  focus  is  not  actually  a point  but 
an  area  of  convergence  of  the  rays.  The  critical  positioning  comes 
about  in  moving  the  receiving  slit  exactly  along  the  radius.  Any 
deviation  will  cause  an  apparent  peak  shift  because  the  angular  relationship 
between  the  direct  beam  (0°20)  and  the  receiving  slit  changes. 

Using  an  automated  receiving  slit  bracket  and  an  automated  residual 

stress  program*replicate  measurements  were  carried  with  and  without 

receiving  slit  movement  on  a 1045-1  sample  which  had  a reasonably  sharp 

211  diffraction  profile.  Using  Cr„  radiation,  15  repeated  laeasurements  were  made 

"or 

using  a three  point  parabolic  fit  to  define  the  diffraction  peak  and  a 
The  automated  residual  stress  program  is  described  in  Chapter  4. 
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statistical  counting  accuracy  of  ± 1.72  MPa  (±  250  psi)  an  average 
value  for  the  stress  of  -164.1  MPa  (-23790  pal)  for  the  parafocusing 
geometry  and  -169.8  MPa  (-24636  psi)  for  the  stationary  slit  method  was 
obtained.  While  the  results  are  remarkably  similar,  the  observed  error 
of  one  standard  deviation  from  the  mean  value  was  * 3.6  MPa  (±  520  psi) 
and  ± 7.4  MPa  (±  1080  psi)  for  the  stationary  slit  and  parafocusing 
geometries  respectively.  The  repositioning  of  the  receiving  slit  was 
very  accurate  and  introduced  only  a small  random  error  of  about  ± 4 MPa 
( ± 600  psi) . The  bracket  uses  a worm  gear  and  dovetail  slide  to  move 
the  receiving  slit  and  is  considered  very  sturdy.  Other,  less  perfect 
systems  such  as  manually  repositioning  the  slit,  are  likely  to  cause 
larger  random  errors  in  the  measurement. 

The  effect  of  systematic  errors  for  the  parallel  beam  method  has 
been  studied  by  Fukura  and  Fujirwara^67^  and  Aoyma.^6®^  The  advantage 
of  this  method  is  that  errors  due  to  sample  displacement  are  minimized 
enabling  simple  sample  alignment.  However,  there  are  several  points  that 
are  not  yet  clear.  The  intensity  provided  by  the  parallel  beam  method  is 
less  than  with  the  focusing  methods  and  instrumental  broadening  is  greater 
which  could  combine  to  decrease  the  precision  of  the  method.  These 
points  are  examined  experimentally  in  Sec.  3.8. 

3.2.3  Instrumental  Errors 

(A)  Sample  Dlsplaceawnt 

The  largest  single  source  of  error  is  displacement  of  the  sample. 

If  the  effective  diffracting  volume  is  not  located  at  the  center  of  the 
diffraction  circle,  there  la  a relative  shift  between  f-00  and  as 


seen  in  Fig.  3.3.  This  has  been  treated  by  Cohen^*^  and  French 

who  derive  a formula  to  calculate  the  peak  shift.  Denoting  AX  as 

the  displacement,  the  equation  for  the  error  in  peak  shift  at  two  tilts 


in  degrees  20  is: 


6(a20)sD  = 


££  AX  cosO  f-i-  . 1 

TT  Lr^  R ' sin ( Off) J 


(3.6) 


where  R ' is  given  in  Eq.  3.3  and  R is  the  goniometer  radius.  The 

P UC 

peak  shift  due  to  sample  displacement  has  been  experimentally  measured 
in  this  study  (see  Sec.  3.7)  verifying  Eq.  3.6  for  the  parafocusing  and 
stationary  techniques. 

Cohen  suggests  a simple  and  precise  method  to  check  for  sample 
displacement  on  a diffractometer.  For  cubic  structures  Cohen  showed: 


*hkl  ao  m _ AX  cos2  6 


rgc  8ln0 


(3.7) 


where  aQ  Is  the  extrapolated  lattice  parameter,  ah^  the  lattice  parameter 

at  a peak  hki,  AX  the  displacement  of  the  sample  off  the  true  center  of 

the  diffractometer  and  R^c  the  redius  of  the  goniometer  circle.  From  a 

plot  of  a^^  vs.  ^ the  slope  equal  to  -a^AX/R^  can  be  obtained  and 

the  displacement  calculated.  A positive  slope  means  the  sample  is 

displaced  too  far  back  (i.e.,  towards  the  back  surface  of  the  sample). 

One  must  only  determine  the  angle  of  three  or  more  Bragg  peaks  and  form 

2 

a plot  of  versus  the  function  cos  0/sin0.  This  method  of  alignment 

is  incorporated  in  the  automated  residual  stress  program  used  in  this 
study;  a sub-program  of  the  package  automatically  determines  the  peak 
positions  and  calculates  the  displacement,  AX,  and  repositions  the  sample. 


FIGURE  3.3  Illustration  of  the  effect  of  sample  displacement  on  the  peak  position  when  the  sample 
is  rotated  t degrees. 
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(B)  Effect  of  t Axis  Mlssettlng 

The  axis  of  inclination  of  the  sample,  f,  must  be  coincident  with  the 
20  axis  or  a peak  shift  will  result  as  shown  in  Fig.  3.4.  Marion 
has  derived  this  error  and  calculates  the  peak  shift  in  degrees  20 
to  be: 

. . ”360  . , sinO  cos 0 (1-cos*) 

6(2°V  -IT-  sln  dH-V)  (3*8) 

where  AX'  is  the  effective  displacement.  This  source  of  error  can  be 

important  on  diffractometers  employing  an  attachment  to  give  the  f 

rotation.  On  the  Picker  diffractometer  the  20  and  u>  axis  rotate  about 

the  same  shaft  so  no  mlssettlng  is  possible.  However,  the  <p  axis  of 

a G.E.  ^-circle  was  also  utilized  occasionally.  The  mlssettlng  of  the 

•k 

(0  axis  with  respect  to  the  20  axis  was  measured  to  be  .05  am  giving  an 
effective  peak  shift  of  .005°20  for  a 45°  tilt  (approximately  3 MPa 
(435  psi)  for  Fe). 

An  estimation  of  the  systematic  error  is  calculated  in  the  automated 
residual  stress  program.  This  includes  the  error  due  to  horizontal 
beam  divergence,  sample  displacement  and  f axis  mlssettlng,  which  are 
the  three  major  errors  (see  Table  3.2). 

3.3  METHOD  OF  PROFILE  LOCATION 

3.3.1  Introduction 

Any  consistent  feature  of  the  line  profile  may  be  utilised  to 
define  a reference  point  from  which  the  position  of  the  profile  is  based. 

*By  mounting  a sample  on  the  ^-circle,  both  the  9 axis  of  the  \ circle  or 
the  tv  axis  of  the  Picker  diffractosMter  could  be  used  to  accomplish  the 
f tilt.  The  mlMsttlng  of  the  £ circle  pea  obtained  by  comparing 
the  peak  position  after  rotating  the  sample  using  each  axis  ror  the  tilt. 


FIGURE  3.4  Illustration  of  the  effect  of  f rotation  axis  not  corresponding  to  the 
20  axis  (point  A is  the  axis  of  20  and  point  0 is  the  axis  of  f). 
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The  position  located  by  any  such  method  need  not  correspond  with 
that  defined  by  any  other  as  long  as  the  reproducibility  is  sufficient 
for  the  measurement.  In  reaidual  stress  measurements  using  X-ray  diffraction, 
the  profile  haa  been  traditionally  defined  in  the  U.S.  by  the  apex 
of  a parabola  fit  to  the  top  region  of  the  curve.  Qualitatively,  this 
gives  equal  weight  to  all  observed  data  points  to  define  the  peak 
location.  It  has  been  suggested^0)  that  the  center  of  gravity  or 
centroid  gives  a more  reproducible  location  of  the  peak,  however,  data 
was  given  only  for  the  422  peak  of  an  annealed  gold  powder  which  is  a 
well  defined  peak  having  a good  peak  to  background  ratio.  In  Japan, 
the  half-value  breadth  and  quarter-value  breadth^^  method  are  the 
adopted  feature  to  represent  the  profile  position.  These  involve  using 
the  midpoint  of  a chord  drawn  through  the  profile  at  a particular 
height.  These  three  methods,  the  peak,  centroid  and  half-value  breadth 
are  compared  for  reproducibility  in  Sec.  3.3.2.  The  peak  is  shown  to 
be  the  most  reproducible  because  the  background  does  not  need  to  be 
determined  since  only  data  near  the  peak  and  not  across  the  whole  profile 
is  usad.  This  implies  that  some  method  must  be  used  to  define  the 
region  of  parabolic  fitting  and  two  such  methods  are  compared  in  Sec.  3.3.3. 
The  formulae  for  a least-squaraa  parabola  is  derived  in  Sec.  3.3.4. 

3.3.2  Comparison  of  the  Precision  of  Various  Measures  of  Profile  Location 
For  the  X-ray  stress  measurement,  the  determination  of  the  diffraction 

angle  by  any  reference  point  on  the  profile  is  aufficlent  because  the 

.*  . b : v.3» '-sr'S  , d , .1  )£/j  ' *rjct  \ . *•/  * « a .)  . rf. 

calculation  of  atrass  is  based  only  on  the  relative  change  of  the  dif- 
fraction angle.  The  precision  of  the  centroid,  half-value  breadth  and 
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peak  positions  ars  compared  using  a position  sensitive  detector  (PSD) 
to  collect  the  date  (see  Sac.  2.5).  The  PSD  system  collects  dats 
across  the  entire  diffraction  profile  slsai ltaneously  so  that  the  data 
for  any  reference  position  can  be  obtained  in  the  same  time  as  for 
any  other  method. 

In  the  half-value  breadth  method  shown  in  Pig.  3.5,  the  diffraction 
angle  of  the  relevant  point  is  determined  by  the  mean  of  two  angles 
given  by: 


So  * *(X-*  ♦ X+*>  <3'9> 

where  is  the  intersection  of  the  profile  curve  and  the  straight 
line  parallel  to  the  background  at  half  of  the  peak  intensity  (excluding 
the  background). 

The  centroid  or  center  of  gravity  of  a diffraction  profile  la 
defined  as:(74) 

Bsucum  and  Aassons^70^  suggest  the  centroid  can  be  located  with  better 
precision  than  the  peak  position  of  the  profile,  although  they  did  not 
make  a direct  comparison.  The  method  of  centroid  calculation  given  by 
Bsucum  and  Ammons  is  based  on  Simpson's  rule  for  a parabolic  approximation 
to  the  area  under  a curve,  deferring  to  Pig.  3.6,  the  centroid,  denoted 

*The  peak  may  be  defined  using  a Cauchy  prefiie,*72*  Osusslan  distribution, (73) 
2nd,  3rd  or  *th  order  polynomials or  a parabola. <2*)  The  parabola 
was  chosen  to  define  the  peak  becauee  it  has  such  general  acceptance  '7^ 
for  stress  measurements. 


X 


FIGURE  3.6  Centroid  calculation 
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X Is: 


where  A. 


and  M. 


k 

I". 

1-1 

k 

K 

i-1 


(3.11) 


X.+h 


- J*  l(x)dx  - 3 [l(X1-h)+4l(X1)+I(X1+h)] 


X1-h 


X.+h 


XI (X) dX  =-  3 [(Xl-h)I(Xi-k)+4X1I(X1)+(Xi+h)T(Xi+h)]  . 


X1-h 


h represents  the  increment  between  data  points  and 


k . (no.  data  points) -3  + ^ 


(3.12) 


since  three  data  points  are  required  for  the  firat  increment  and  two 
for  the  rest. 

Comparisons  of  the  reproducibility  of  the  least-squares  parabola 
(see  Sec.  3.3.4  for  the  formula),  half-value  breadth  and  centroid  were 
made  on  the  1090-1  and  1045-2  samples  using  the  CrKa  211  peak.  The 
automated  residual  stress  program  described  later  for  the  PSD  was 
modified  to  Include  the  calculations  required  in  Eq.  3.9  for  the  half- 
value  breadth  and  Eq.3.11  for  the  centroid.  A linear  background  was 
determined  by  fitting  a straight  line  through  50  data  points,  25  on 
each  aide  of  the  tails  of  the  diffraction  profile.  The  intensity 
corrections  (Lorents-polarizatlon  and  absorption)  were  laade  after  the 
background  had  bean  subtracted  and  prior  to  determining  the  centroid 
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and  position  of  the  half-value  breadth.  The  least-squares  parabola  was 
fit  to  the  top  15  pet  of  the  profile  without  first  subtracting  the 
background. 

The  results  of  the  ten  replicate  measurements  are  given  in  Table 
3.3.  The  two-theta  position  represents  the  average  calculated  value 
over  the  ten  measurements  with  one  standard  deviation  from  the  mean 
being  given  in  parenthesis.  For  the  1090-1  sample,  having  a sharp 
diffraction  profile  and  low  background,  the  centroid  actually  has  a 
smaller  variance  than  either  of  the  other  two  methods.  For  the  1045-2 
sample  which  has  a poor  peak  to  background  ratio,  the  parabola  has  the 
lowest  variance.  For  both  samples,  the  half-value  breadth  gave  the 
poorest  reproducibility. 

In  X-ray  residual  stress  measurements  the  peak  to  background  ratio 
is  usually  poor,  especially  for  hardened  materials.  In  the  centroid  and 
half-value  breadth  reference  points,  the  value  Is  obtained  from  a 
truncated  portion  of  the  profile  determined  by  the  background  subtraction. 
Thus,  the  background  subtraction  introduces  an  error  (often  referred 
to  as  the  truncation  error'  ').  For  the  least-squares  parabola,  even 
if  background  is  subtracted  in  determining  the  region  of  curve  fitting 
(necessary  only  when  a few  data  points  are  obtained,  see  Sec.  3.6.3) 
it  does  not  affect  the  precision  because  the  entire  profile  is  not  used 
minimizing  the  truncation  error. 

3.3.3  Region  of  Parabolic  Fit 

The  region  in  the  vicinity  of  the  peak  that  can  be  used  for  the 
parabolie  fit  is  not  wall  defined.  If  it  la  assuamd  the  true  profile 
can  be  fit  by  a quasl-Lorentztlan  function  (which  is  now  becoadng  quite 
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TABLE  3.3 

PRECISION  OF  VARIOUS  MEASURES  OF  PROFILE  POSITION 
(10  measurements) 


Time 

(sec) 

FWHM* 

(°20) 

Half-Value  Breadth 
(°  20) 

Centroid 

(°20) 

Parabola 

(°20) 

50 

.45 

156.149  (±  .021)** 

156.096  (±  .011) 

156.186  (±  .016) 

100 

3.45 

155.336  (±  .064) 

155.396  (±  .085) 

155.413  (±  .020) 

* 

FWHM  is  the  full  width  at  half  of  the  maximum  Intensity. 

** 

The  term  in  ( ) represents  one  standard  deviation  from  the  average 
position  over  the  10  measurements. 
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popular  In  profile  fitting  powder  patterns  for  structure  analysis  ) , 
Thomsen  and  Yap^62^  have  shown  that  a parabola  Is  valid  for  V less  than 
.32  where: 


V 


2<20p  - 

W 


20min> 


(3.12) 


Here  W is  the  full  width  at  half  maximum  height  of  the  diffraction  pro- 
file, the  peak  position  is  20^  and  the  minimum  value  of  20  lying  on 
the  parabolic  curve  is  20^^.  While  being  theoretically  justified 
Eq.  3.12  still  depends  on  the  profile  being  symmetric  about  the  peak. 

Koistinen  and  Marburger ^26^  suggest  that  data  points  for  the 
parabolic  fitting  should  be  chosen  which  have  intensities  at  least  857. 
of  the  maximum  intensity.  This  rule  has  gained  wide  acceptance  in  X-ray 
residual  stress  analysis  although  it  is  purely  empirical. 

The' top  15  pet'  rule  is  commonly  used  in  X-ray  residual  stress 
analysis  because  it  is  simple  to  apply.  Only  data  near  the  top  region 
of  the  peak  need  be  recorded  and  processed.  The  formula  given  by 
Thonsen  and  Yap  requires  knowledge  of  the  breadth  of  the  peak  at  half 
maximum  intensity  which  demands  accumulating  data  across  the  entire 
profile  and  determining  the  background. 

The  two  rules  for  Judging  the  region  of  the  parabolic  fit  have 
been  applied  to  experimental  data.  Table  3.4  tabulates  the  peak 
determined  from  a three  point  parabolic  fit  in  column  4 and  the  minimum 
20  setting  to  be  on  a parabola,  as  calculated  by  each  method  in  the 
next  two  columns.  Column  2 is  the  full  width  at  half  the  maximum 
intensity  while  column  3 records  the  ratio  of  the  peak  intensity  to 
the  background  intensity. 
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For  the  first  two  samples  listed,  both  having  reasonably  sharp 
diffraction  profiles,  the  'top  15  pet'  rule  of  Koistinen  and  Marburger 
corresponds  well  with  that  of  Eq.  3.12.  This  similarity  breaks  down 
as  the  profile  becomes  broader,  as  depicted  by  the  remaining  four  samples. 
For  these  samples  exhibiting  increasingly  broader  diffraction  profiles, 
the  decreasing  similarity  between  the  'top  15  pet*  rule  and  Eq.  3.12 
corresponds  to  a decreasing  peak  to  background  ratio  (column  3) . The 
background  was  not  subtracted  in  determining  the  top  15  pet  region 
resulting  in  the  region  being  substantially  larger.  If  a linear  back- 
ground is  subtracted  before  applying  the  '15  pet'  rule,  the  regions 
correspond  quite  closely,  as  seen  in  column  7.  This  background  correction 
need  not  be  exact.  A quick  estimate  of  the  background  obtained  at  one 
20  position  away  from  the  peak  is  sufficient  since  this  is  only  concerned 

with  determining  the  region  of  curve  fitting.  This  procedure  is  incorporated 
in  the  computer  program  STRESS  described  in  Chapter  4. 

3.3.4  Derivation  of  Peak  Location  using  Parabolic  Fit 

It  is  assumed  that  the  data  is  obtained  at  an  odd  number  of 
observation  points  and  is  also  taken  in  equal  increments  of  20. * The 
data  can  then  be  said  to  be  measured  at  2n+l  points,  the  center  point 
being  taken  as  a working  origin.  The  data  accumulated  at  each  jth  point 
is  Ij,  the  power  in  counts  per  second.  If  the  20  Increment  between 
data  points  is  designated  as  6,  the  parabola  is  defined  by: 

a + bfij  + c«2J2  - Ij  (3.13) 


★ 

Data  taken  using  the  PSD  is  in  equal  Increments  of  distance  which  does 
not  translate  exactly  into  equal  20  increments.  The  error  in  using  a 
linear  calibration  constant  to  transfer  distance  along  the  PSD  into  °20 
is  minimal  for  rasidual  strass  swasurements  (see  Sec.  2.5.5). 


39 


and  the  least-squares  solution  is  found  by  minimizing: 

n 

S = ^ (a  + bftj  + C62j2  - Ij)2  (3.14) 

J—n 

with  respect  to  the  unknown  parameters  a,  b,  and  c.  To  solve  Eq.  3.14, 
the  partial  derivatives  of  S with  respect  to  each  parameter  are  set 
equal  to  zero.  This  gives  three  equations  for  estimating  the  values  of 
a,  b,  and  c to  minimize  S. 

For  example: 


* 0 - 2 ^ (*  + bfij  + cfiZj2  - 1^) 

J—n 


(3.15) 


• 1 J°  + b8  1 J + c«2  'L  J2 

j—n  J—n  J—n  j—n 
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where  n * I1 
i U J 

J—n 

referred  to  as  the 
Continuing  in 


n 

and  * 6 1+1  ^ J*Ij  . The  term  will  be 

J—n 

i moment  of  the  observed  line  profile, 
this  manner  the  set  of  equations  are: 


2 -1 

n a + 0 + n,6  c - 6 M ■ 0 

oz  o 

0 + t^ftb  + 0 - 6'2Ml  - 0 

2 -3 

n2a  + 0 + n^6  c - 6 Mj  - 0 
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Wilson^63^  has  derived  the  statistical  variance  for  a least-squares 

parabola  fitted  to  data  on  X-ray  peaks  collected  for  fixed  time  at 

each  observation  point.  Wilson's  work  is  extended  here  to  include 

both  fixed  time  and  fixed  count  methods  of  data  collection  to  obtain 

the  counting  statistical  error  in  a peak  location.  Propagation  of  this 

error  to  the  error  in  the  calculated  stress  is  presented  for  both  the 
2 

'sin  f ' and  'two  tilt'  methods  of  X-ray  residual  stress  measurement  in 
Sec.  3.5. 


3.4.2  Derivation  of  Statistical  Variance  in  the  Peak  Location 

The  variance  of  a general  function,  X,  of  several  variables, 
(76) 


X— f (x^ , X£ , • • • , x^)  is: 


,.2 


V<X>  ■ Csp  v<xl>  * (sp  V<V+-+<E7  i^C0V<V2>+- 

(3.18) 

+ higher  order  differential  terms.  Here  V(X)  is  the  variance  of  X and 
COV(x^ ,Xj)  is  the  covariance  between  variables  x^  and  x^.  Neglecting 
the  higher  differentials  this  can  be  rewritten  as: 


(3.19) 


v<x)  " 

If  the  function  is  linear,  X ■ V a^x^,  the  standard  error  is  given  by 


■? 


V (X)  - ) a1V(x1) 


(3.20) 
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Since  the  count  rate,  1^  s is  « random  variable  having  a finite 
variance,  all  quantities  obtained  by  manipulating  it  also  have  a finite 
variance.  Since  the  measured  Intensities  are  statistically  independent: 


COV(xt>Xj)  * 0 for  i i J 

COV(x1,Xj>  - V(xt)  for  i - J 


(3.21) 


Eq.  3.19  then  becomes: 


v2 

V(X)  V(xj>  * <3*22> 

J ' 

This  equation  is  reasonably  accurate  for  any  nonlinear  function  providing 
the  standard  deviation  of  each  variate  is  small,  say  20 X of  the  mean.^?6) 

Applying  Eq.  3.22  to  Eq.  3.17  and  substituting  for  x^  the  counting 
rate,  I ^ : 


(3.23) 


The  variance  of  the  peak  is  dependent  on  the  variance  of  the  count  rate, 

Ij,  which  depends  on  how  the  data  Is  accumulated.  The  two  cotason 

methods  of  determining  the  counting  rate  at  the  diffractometer  setting 

2®.  the  fixed  time  and  fixed  count  techniques.  Defining  t.  as  the 

J J 

time  of  data  accumulation  and  as  the  accumulated  counts: 


V' 

for 

fixed  time 

c/tj 

for 

fixed  counts 

(3.24) 


A3 


Coafclnlng  Eq.  3.23,  3.25  and  3.29  will  glva  the  error  in  peek 
location  due  to  counting  atatlatlca  alone. 

For  accuMilatlon  of  dace  by  fixed  tiaw  (at  each  obaervation  point 
j)  the  foraula  ia; 

I ces.*j 
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° <2V 


4.  2 ,2  n 

o ln,-n  n ) r 3 2 2 2 3 2 

2.  .2,  ..  I l <V\1  +(V  VW‘  J‘"2«  V 1J/t 


4W  VW  j--n 


(3.30) 

where  o(20^)  represents  the  standard  deviation  in  the  peak  position. 

The  assumption  of  a parabolic  profile  implies  symmetry  about  the  peak 
position.  The  term  represents  an  I1*1  moment  of  Intensity  about 
the  peak.  Wilson suggests  that  because  the  profile  has  been  con- 
sidered syranetrlc  (at  least  in  the  region  fit  by  a parabola)  the  odd 
moments  of  Intensity  will  be  of  less  liaportance  (for  a perfect  para- 
cola  the  odd  moments  are  zero)  reducing  Eq.  3.30  to: 


a (20  ) - 


4 2 .2  n 

6 (n2-nn4)  ..  2..  v 2.4.2 


4‘t-n2(n26  M^n^) 


4 ).  <”2«  VW  5 J >) 


,4,  2 \ 2 

6 (n2-non4) 


2 2 

4*t*n2(n2C 


.4  \ 


1 * 6 i i~l‘ 

J*-n 


.8/  2 v2  „ 

6 (n2-nQn4) 

2 — r 2 

4-t»n2(n26  l^-n^) 


(3.31) 


In  the  same  manner,  for  data  accumulated  by  fixed  counts  at  each 


-2(2V ' ~r— j 


.4.  2 %2 

« (n2-non4) 


4n2(n25  M^nMj)  J--n 


(3.23) 


Neglecting  the  odd  momenta  of  Intensity,  Eq.  3.32  reduces  to: 


=2<2y 


,8,  2 2 n 

6 (n2-non4)  2_  2 

“T  2 / J 

4CW  VW  i-n  1 


(3.33) 


3.4.3  Experimental  Confirmation  of  the  Statistical  Formulae  for  Peak 
Location 

The  counting  statistical  errors  for  the  peak  location  were  examined 
with  the  211  peak  using  three  samples  to  cover  a range  of  breadths;  the 
1090-1,  1045-1,  and  TBA  G-5  specimens. 

Three  observation  points  were  chosen  from  a chart  scan  of  the  profile 
and  were  within  the  top  15  pet  of  the  peak  intensity.  Data  was  repeatedly 
accumulated  using  fixed  counts  at  these  same  observation  points  15  times 
for  each  sample  and  the  peak  position  and  statistical  error  (using  Eq.  3.32) 
were  calculated  employing  a computerized  system.*  The  results  are 
given  in  Table  3.5.  The  peak  position  is  the  average  value  of  the  15 
measurements.  The  standard  deviation  from  the  mean  peak  location, 

S(29p)  is  given  in  column  3.  When  compared  with  the  average  counting 
statistical  error  (one  standard  deviation)  given  in  column  4,  it  is 
seen  that  the  correspondence  is  remarkable  thereby  showing  that  Eq.  3.32 
for  fixed  count  analysis  predicts  the  proper  counting  error.  The  slight 
difference  is  probably  due  to  mlssetting  of  20  because  the  limit  of 
gearing  accuracy  on  the  Picker  diffractometer  is  .002°20. 

★ 

The  program  as  described  in  Chapter  4 was  modified  slightly  to  use 
the  same  angular  position  for  tha  data  points  during  the  replicate 
tests.  As  will  be  seen  in  Chapter  4,  the  program  normally  selects  these 
positions  and  in  each  maasurement  the  different  positions  would  affect 
these  tests. 
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TABLE  3.5 

ERROR  ANALYSIS  FOR  REPEATED  PEAK  DETERMINATIONS  Cr  - 211  PEAK 

(15  Measurements) 


Sample 

S(20  ) 
°2ok 

o(20p) 

°20P 

1090-1 

156.184 

.0022 

.0021 

1045-1 

155.916 

.0023 

.0021 

TBA  G-5 

155.282 

.0112 

.010 

*(20  ) repretent*  the  mean  value  of  the  peak  poaitlon. 
p 

S(20  ) rapreaenta  the  observed  atandard  deviation  froai  the 
p mean  value. 

o(20^)  ia  the  counting  atatiatical  error. 
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3.4.4  Experimental  Confirmation  of  the  Approximate  Formulae  for  the 
Variance  in  Peak  Location 

By  asaumlng  the  actual  data  to  be  symmetric  about  the  working  origin, 
i.e.,  the  center  of  the  curve  fitting  region,  Eq.  3.31  (for  fixed  time) 
and  Eq.  3.33  (for  fixed  counts)  were  derived,  A test  of  these  appro- 
ximations were  made  by  analyzing  data  obtained  from  the  Cr„  211  / 

diffraction  peak  on  steel  samples.  The  tests  were  carried  out  for  both 
and  ^-45°  (after  tilting  the  peak  may  be  more  asymmetric)  using 
stationary  slit  geometry.  The  1090-1  and  TBA  G-5  samples  ye're  used, 
the  first  exhibiting  excellent  separation  with  the  second  having 

no  separation.  Table  3.6  reproduces  the  results  for  the  fixed  count 
method.  The  nuu&er  of  data  points  is  given  in  Column  3 and  the  step 

increment  between  points,  6,  is  given  in  column  4.  The  breadth  at  half 

/ 

the  maximum  Intensity  is  given  in  column  5. 

Three  data  points  from  a peak  of  any  shape  can  form  a perfect 
parabola.  Therefore,  one  would  not  expect  any  difference  in  the  results 
between  the  exact  formula  and  the  approximation  (based  on  perfect 
aynnetry  of  the  actual  data)  using  three  data  points.  This  is  seen  to 

be  so  in  Table  3.6  for  both  samples.  With  a multi-point  fit,  any 

/ 

/ 

asymmetry  will,  be  more  evident.  From  Table  3.6,  the  difference  between 
the  results^  calculated  by  both  formulae  is  negligible  for  the  multipoint 
fit  implying  that  the  actual  profile  is  aymmatric  about  the  peak 


location  after  the  data  has  been  corrected  for  the  20-dependent  intensity 

factors.  The  exact  formula  is  used  in  the  computer  program,  STRESS  described 
in  Chapter  4. 

3.3  STATISTICAL  VARIANCE  IN  STRESS 

The  variance  in  the  measured  peak  position  results  in  an  error  in 


J 
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the  calculated  residual  stress.  This  error  will  be  derived  for  both 
2 

the  'sin  method  and  the  'two  tilt'  method. 

3.5.1  *Sln2f’  Technique 
2 

In  the  'sin  s'  technique  the  residual  stress,  a , is  determined 

2 

from  the  slope  of  the  interplanar  spacing  d versus  sin  *.  The 

5P»f 

variance  in  d due  to  random  errors  in  the  measured  intensities  will 
cp»  f 

introduce  an  error  in  the  surface  stress.  The  variance  in  d must  be 

<P  .* 

determined  from  the  variance  in  20  . 

P 

Writing  Bragg* 8 law; 


d - X/2sin0 


(3.34) 


Applying  Eq.  3.20  to  the  relation  in  Eq.  3.34: 


o2(d) 


o2(0) 


' £3?*  ^ 


(3.35) 


where  20  is  in  degrees  and  a (20)  is  given  by  any  one  of  Eqs.  3.30  - 3.33. 

The  surface  stress  is  given  by  m',  the  regression  coefficient  of 
2 

d ^ vs.  sin  f.  The  regression  coefficient  of  a linear  least  squares 
line  la  given  by:^"^ 


>•  - 


E(*r*) (y.-y) 

1 1 

— _ .■ 

Z(xt-x) 

.»(.f  r> 


(3.36) 
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where  is  the  independent  veriable  and  y^  is  the  dependent  variable. 

2 

The  bar  denotes  averages,  in  terms  of  d . and  sin  i , the  regression 

«P»f 

coefficient  is  expressed  as: 


where 


ZC.in  f-ait?^i)  (d  . “d  . ) 


- L 


cp. t tp.r 


Z(ain2t-sT 

♦ 

sin2*  - Z sin%/N 

♦ 


(3.37) 


«L  . - E d /N 
*P»?  . <P»T 


and  N is  the  total  number  of  observations. 

The  value  of  * also  has  finite  error  due  to  alignment  and  mechanical 

2 

motion.  The  function  sin  * varies  slowly  so  small  errors  in  * have 

2 

negligible  effect  on  sin  *.  Thr  expression  for  the  regression  co- 
efficient, Eq.  3.37  is  then  a lineax  function  with  the  dependent  vari- 
able having  a predictable  error.  By  application  of  Eq.  3.20,  the  variance 
of  m*  Is; 


Z(eln2t-«in^)2  . „2(d  ) 

2 * f 

o (•’)  • — — — 


rZ(aln2t-ain%)2"|2 

t 


(3.38) 


where  g (d  ) is  given  by  Eq.  3.35.  Through  relation  3.20,  the  variance 
v>»? 

in  the  stress,  g , (given  by  Eq.  1.5),  is: 

V 


2 

o2(0j  - a*  (*')/&, 


(3.39) 


where  g (m*)  la  given  by  Eq.  3.38. 
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3.5.2  'Two  Tilt'  Technique 

Application  of  Eq.  3.19  to  the  'two  tilt'  method  (given  by  Eq.  1.6)  yields 


r.W<-#>]  (3.40) 

o T 

where  the  stress  constant  K is  given  by: 


(3.41) 


If  the  peak  shift  is  small,  the  two  tilt  equation  can  be  expressed 
directly  in  terms  of  20  by  Eq.  1.7.  The  variance  in  the  surface  residual 
stress  is  given  by: 


* K,202( 20oHo2(26^)]  (3.42) 

where  the  variance  in  peak  position  is  given  by  Eq.  3.30  or  Eq.  3.32. 
Small  errors  In  20  do  not  introduce  significant  error  in  the  stress 
constant  K* , given  by  Eq.  1.8. 

3.5.3  Comparison  of  Statistical  Formulae  for  "IVo  Tilt'  Method 
Kelly  and  Short^®^  derived  an  equation  for  the  error  in  the 
residual  stress  as  determined  using  the  'two  tilt'  method  due  to 
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counting  statistics.  Using  Eq.  1.8  and  writing  the  peak  location 

directly  in  terms  of  the  counts  at  each  of  three  data  points  (therefore 

using  a three  point  fit)  Kelly  and  Short  obtained  an  equation  relating 

the  stress,  a , directly  to  the  intensity  obtained  at  the  three  obser- 
«P 

vation  points.  They  applied  Eq.  3.22  to  this  formulae  to  calculate 
the  counting  error.  Defining  (1-P)  as  the  relative  difference  in 
intensity  between  the  peak  position  and  the  outer  observation  points, 
and  N as  the  total  counts  at  each  point,  they  obtained,  for  fixed 
count  data  accumulation: 

* 

o(ow)  - K-A20.yw  (2(1-P)/N)  (3.43) 

where  K is  the  stress  constant  in  Eq.  1.8.  This  equation  is  valid  only 
for  a three  point  parabola  and  the  'two  tilt'  method  of  residual  stress 
analysis.  Their  derivation  assumes  the  angular  increment  between  data 
points  to  be  identical  for  both  f inclinations  and  also  the  same  total 
counts  are  achieved  at  both  f tilts. 

Table  3.6  compares  results  obtained  from  analyzing  the  expected 
counting  error  with  Eq.  3.42  and  Eq.  3.43  from  data  taken  on  three  steel 
specimens.  The  data  was  taken  in  the  top  157.  region  using  fixed  counts. 
The  samples  were  chosen  to  cover  a range  of  peak  breadths  and  residual 
stress  levels. 

The  counting  errors  calculated  by  each  aquation  are  given  in  Table 
3.7  and  are  similar,  although  exact  agreement  is  not  expected  because 
of  the  assumptions  made  by  Kelly  and  Short.  The  results  do  indicate 
the  statistical  error  from  both  methods  are  comparable,  however,  the 
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TABLE  3.7 

COMPARISON  OF  THE  STATISTICAL  ERROR  IN  RESIDUAL  STRESS 
AS  CALCULATED  IN  THIS  STUDY  AND  BY  KELLY  AND  SHORT<78) 

Cr  radiation  - 211  peak  - three  point  parabolic  fit 
Ko 


Sample 

* 

6 

°20 

B 

cr(cffl),MPa 
Eq.  3.43 

o(o8)  ,MPa 
Eq.  3.42 

1090-1 

+12.8  (+1860  psl) 

.11 

.155 

±1.97(140  pal) 

±.875  (127  psl) 

1045-1 

-176  (-25,600  psl) 

.295 

.159 

±3.49  (506  pel) 

±3.14  (456  psl) 

1045-2 

-400  (-57,950  psl) 

2.11 

.172 

±8.17  (1185  psl) 

±8.96  (1300  psl) 

6 la  the  step  increment  between  the  angular  positions  of  the  data  points. 
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lack  of  generality  of  the  Kelly  and  Short  analyaia  limits  its  application. 

Their  derivation  is  based  on  the  two  tilt  method  for  the  residual 

stress  determination  and  assumes  the  peak  shift  (°20)  and  stress  are 

proportional.  This  is  valid  only  for  a small  peak  shift.  The  statistical 

analysis  as  developed  in  Sec.  3.4.2  enables  the  statistical  error  to  be 

calculated  for  either  the  exact  or  approximate  'two  tilt'  formulae 

2 

(Eq.  3.40  and  Eq.  3.42,  respectively)  and  for  the  'sin  ' method  (Eq. 

3.39). 

3.5.4  Automated  Residual  Stress  Program 

A complete  computer  controlled  residual  stress  analysis  program 
was  developed  during  this  study.  An  in  depth  discussion  of  the  program 
is  given  in  Chapter  4 , but  because  the  program  is  used  in  the  remaining 
sections  of  Chapter  3 the  program  is  described  briefly  here  to  facilitate 
the  reader's  understanding. 

The  residual  stress  program  incorporates  the  following  features 
which  are  choser.  by  the  user  by  means  of  an  initial  dialog: 

1)  The  experimental  parameters  are  input  (approximate  peak 
location,  divergent  slit,  sample  displacement  and  f-axls  missettlng) 
to  calculate  geometric  aberrations. 

2 

2)  Either  the  'two  tilt'  or  'sin  methods  of  residual  stress 
analysis  may  be  applied,  with  a normal  detector  or  the  PSD. 

3)  The  Marlon-Cohen  technique  is  Implemented  automatically  if  I 

2 

a least-squares  fit  to  d vs.  sin  f Indicates  oscillations. 

4)  The  operator  apecifies  the  desired  accuracy  in  tanas  of 
degrees  20  for  each  peak  or  in  stress  and  the  counting  strategy  is 
determined  taking  into  account  both  statistical  and  geometric  errors. 
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5)  The  peak  is  fit  to  a parabola  by  a least-squares  procedure. 

The  operator  specifies  the  number  of  data  points  to  be  fit. 

6)  If  desired,  background  is  measured  and  subtracted. 

7)  Either  parafocusing  or  stationary- slit  geometry  can  be  chosen. 

8)  An  estimate  of  the  overall  time  of  analysis  is  given  after 

a preliminary  scan  of  the  peak.  The  operator  can  then  change  the  error 
if  the  time  is  too  long. 

Five  successive  step  scans  are  involved  in  the  procedure  for 
determining  a peak  position,  as  described  below.  The  data  is  collected 
in  fixed  count  mode  and  corrected  for  deadtime  and  angular  dependent 
intensity  factors. 

1)  A step  scan  is  made  in  large  increments  (~  ,2°20)  set  by  the 
operator,  accumulating  1000  counts  at  each  position  from  the  initial 
20,  until  a count  rate  of  less  than  90  pet  of  the  maximum  is  obtained. 

2)  TWo  step  scans  follow  in  smaller  increments  (set  initially  by 
the  user  and  typically  0.05°)  down  each  side  of  the  peak  to  locate  two 
angles,  20^  and  20^  at  85  pet  of  the  maximum  intensity.  (The  background 
measured  at  an  angular  position  specified  by  the  user,  may  be  subtracted 
from  the  data  in  determining  the  two  end  points.  This  procedure  serves 
to  define  the  region  of  parabolic  fit  more  accurately.) 

3)  A third  step  scan  at  the  two  angles  from  Step  2 and  the 
central  angle  between  them  is  made  with  a preset  count  of  5000.  A 
three  point  parabola  is  fit  to  the  data.  Then  angles  are  calculated 
for  the  desired  number  of  data  points  such  that  the  central  data  point 
will  he  very  close  to  the  center  of  the  parabolic  fit.  This  minimizes 
odd  order  terms  in  the  error  equations.  These  preliminary  scans  typically 
take  60  seconds • 
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A)  A step-scan  at  the  final  angular  positions  for  a preset  1000 
counts  is  used  to  back  calculate  the  necessary  counts  from  Equation  3.32 
to  obtain  the  desired  precision.  Steps  3 and  A constitute  a multiple 
pass  procedure  which  serves  to  improve  the  reliability  of  peak  location 
and  minimize  the  time  required  for  a given  precision. 

5)  Hie  final  data  is  acquired  for  the  calculated  preset  counts 
at  each  angular  position. 

3.6  EXPERIMENTAL  DETERMINATION  OF  THE  PRECISION  OF  RESIDUAL  STRESS  ANALYSIS 
3.6.1  Introduction 

A complete  statistical  analysis  of  the  peak  location  using  a least- 
squares  parabolic  fit  to  the  top  of  a diffraction  profile  has  been 
presented.  While  useful  in  determining  random  errors  due  to  counting 
statistics,  other  factors  also  contribute  to  the  precision  of  the 

o 

measurement.  The  experimental  technique  (i.e.,  'two  tilt'  vs.  'sin  f' 

methods)  and  accuracy  of  the  equipment  influence  the  overall  precision. 

The  actual  precision  has  been  examined  in  this  atudy  under  varying 

conditions.  The  reproducibility  of  residual  strass  determination  using 

a least  squares  parabolic  curve  for  peak  location  is  compared  to  that 

using  a three  point  parabola.  The  precision  of  the  'two  tilt'  method 
2 

and  'sin  f'  method  are  examined  and  compared  to  each  other.  These 
treatments  on  precision  have  apparently  never  before  been  carried  out. 

The  reproducibility  of  residual  stress  determination  using  a 
position  sensitive  detector  (PSD)  for  data  accumulation  is  alao  examined. 


The  precision  is  dependent  on  the  resolution  and  stability  of  the  PSD 
and  these  are  shown  to  be  quite  adequate  for  residual  stress  analysis 
by  X-ray  diffraction. 
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The  precision  of  the  residual  stress  measurement  Is  dependent  of 
the  method  of  paak  location.  Traditionally  in  X-ray  stress  steasurements 
using  a diffractometer,  the  apex  of  a three  point  parabola  fit  to 
the  top  region  of  the  diffraction  profile  is  used  to  define  the  Bragg 
angle  as  has  been  mentioned  several  times,  it  was  shown  in  Sec.  3.4.3 
that  the  region  near  the  maximum  of  the  profile  is  symmetric  about  the 
peak  position  after  Lorentz-polarlzatlon  and  absorption  factors  are  made, 
verifying  the  practicality  of  the  three  point  parabola  compared  to  a 
least-squares  approach.  Jatczak  and  Boehm  v conclude  that  a sore 
elaborate  fitting  of  2nd,  3rd  and  4th  order  curves  to  five  data  points 
produces  good  agreement  with  the  three  point  parabola  method  and  hence 
was  not  worth  the  extra  time  needed  in  accumulating  the  data.  However, 
Marion*30^  found  a least-squares  parabolic  fit  to  10  to  20  data  points 
gave  better  reproducibility  than  Just  the  three  point  fit.  This 
apparent  contradiction  is  important  to  resolve  because  the  curve  fitting 
procedure  will  affect  the  precision.  It  is  obvious  that  many  data  points 
are  preferable  to  only  a few  given  unlimited  time.  For  a fixed  total 
time,  however,  this  is  not  necessarily  true.  Yap*67^  concluded  that  if 
the  profile  is  a perfect  parabola  the  optimum  procedure  for  collecting 
data  is  to  spend  most  of  the  time  at  the  two  end  points  of  the  parabola 
and  only  a short  time  at  the  peak  in  a three-point  fit.  Since  random 
errors  dictate  that  the  observed  data  la  not  a perfect  parabola  the  effect 
of  the  number  of  observation  points  on  the  precision  was  experimentally 
f®*t*d  in  this  study.  The  total  time  for  data  accumulation  was  fixed 
and  tests  were  run  using  the  'two  tilt’  method  in  the  following  ma 


inner: 
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1)  Fix  che  total  time  spent  accumulating  the  data. 

2)  Determine  the  time  at  each  point.  Counting  can  be  accomplished 
in  either  fixed  time  or  fixed  count  mode. 

3)  Repeat  the  measurement  until  enough  results  have  been 
accumulated  to  check  the  reproducibility. 

4)  Duplicate  steps  2 and  3 for  a different  number  of  data  points, 
keeping  the  total  time  of  measurement  constant. 

By  varying  the  number  of  data  points,  the  effect  of  this  number  on 
the  precision  can  be  assessed. 

3.6.2  Experimental  Procedure 

An  automated  diffractometer  system  was  used  to  accomplish  the 
repeated  measurements.  The  system  is  comprised  of  a Picker  diffractometer 
and  scintillation  detector  interfaced  to  a PDP-8E  computer.  A position 
sensitive  detector  was  also  used  in  conjunction  with  the  diffractometer. 
The  features  of  the  software  programming  are  described  in  Sec.  3.5.4. 

The  only  modification  of  the  program  for  these  tests  concerned  the  pro- 
cedure for  determining  the  fixed  counts  accumulated  at  each  observation 
point.  In  its  usual  form,  the  computer  program  allows  the  operator  to 
chose  the  desired  number  of  data  points  and  determines  the  necessary 
counting  time  to  achieve  a designated  atatistical  precision.  Since  the 
total  time  of  data  accumulation  must  be  fixed  for  each  peak  determination 
in  this  analysis,  the  following  procedure  was  substituted.  Defining 
the  total  time  as  T,  the  fixed  counts  at  aach  observation  point,  J, 
as  CT,  and  the  power  in  counts/aacond  at  each  j as  IN(J),  tha  avarage 
power  tiswa  total  time  will  yiald  the  total  counts.  Dividing  this  by 
the  number  of  data  points,  N,  yields  an  average  fixed  counts.  Therefore: 


where  IN(J)  represents  en  average  Intensity  In  CPS.  The  individual  IN(j) 
are  determined  initially  from  a rapid  step  scan  used  to  determine  the 


position  of  the  observation  points  through  the  use  of  the  'top  15  pet 


rule.  The  total  time,  T,  was  input  and  the  program  determined  CT 


For  the  position  sensitive  detector,  the  total  time  is  fixed  by 
the  multichannel  analyzer  which  acts  as  a peripheral  storage  device  for 
the  data  so  it  was  not  necessary  to  alter  the  usual  computer  program. 
The  samples  were  chosen  to  cover  a broad  spectrum  of  residual 


and  breadths  of  diffraction  profiles 


3.6.3  The  Effect  of  the  Number  of  Data  Points 


Ten  repeated  stress  measurements  were  made  on  each  of  the  steel 


samples  for  each  number  of  observation  points.  Three, 

five,  seven,  eleven  and  fifteen  observations  in  the  region  of  curve 
fitting  were  used.  The  diffractometer  allgnawnt  was  checked  and  the 
sample  positioned  for  minimum  sample  displacement  prior  to  beginning  the 


measurements  on  each  sample.  No  instrumental  changes  were  made  during 
the  replicate  testing. 

The  'two  tilt'  method  using  f angles  of  (fi  and  45°  was  employed 
with  a stress  constant  of  593  MPa/°2Q  (86000  psl/°20).(7)  The  stationary 
slit  or  non-focusing  technique  was  used  so  as  not  to  include  random 
errors  due  to  repositioning  Of  the  receiving  slit  (see  Sec.  3.2.2  for 
more  details  on  the  effect  of  slit  positioning). 


Table  3.8  sumnarlses  Che  effect  of  the  number  of  data  points  employed 
to  determine  the  peak  maximum  on  the  precision  with  a normal  detector 
and  the  two  tilt  procedure.  The  number  of  data  points,  column  2 and 
the  total  time  of  data  accuaailatlon  for  the  two  tilts,*  column  3,  were 
set  In  the  initial  dialog  with  the  computer  program.  The  residual  stress 
given  In  column  4 is  the  average  value  over  10  measurements.  The  observed 
error,  one  standard  deviation  from  the  mean  stress  for  the  10  replicate 
measurements  is  given  in  column  5 and  the  statistical  counting  error 
for  one  standard  deviation  (Eq.  3.42)  is  given  in  the  last  column. 

The  data  for  the  1090-1  sample  which  exhibits  a sharp  211  diffraction 
profile  with  excellent  K^-K^  separation  illustrates  that  for  such  a 
profile,  the  number  of  data  points  has  little  effect  on  the  precision. 

The  counting  error  is  slightly  less  than  the  observed  precision,  as 
expected.  For  this  sample  the  precision  is  approximately  ±1.4  MPa 
(±200  psi),  better  than  that  usually  reported  in  the  literature  (±  10  MPa)*7) 
because  the  peak  was  very  sharp. 

As  the  breadth  of  the  diffraction  profile  increases  the  peak 
position  becomes  more  sensitive  to  small  fluctuations  In  the  data  and 
the  counting  statistics  for  the  three  point  fit  do  not  predict  the  true 
err®r»  But  as  the  number  of  data  points  In  the  least-squares  fit 
Increase,  the  observed  precision  once  again  approaches  the  counting 
error.  These  results  show  that  counting  statistics  are  not  a good 
measure  of  the  precision  when  using  only  three  points  for  a parabolic 
fit  on  samples  having  broad  profiles.  Only  where  ttm  number  of  data 
points  Is  Increased  Is  the  observed  precision  represented  by  this 

*"  ———————— 

The  total  time  Is  for  the  final  step  scan  only  and  does  not  Include 
the  time  spent  in  the  preliminary  scans  as  this  Is  Independent  of  the 
number  of  data  points. 
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TABLE  3.8 

RESULTS  OF  FIXED  TI*  TESTS 
(TOO  TILT  METHOD) 


Number  of  Total  Time 
Data  Points  (sec) 


Stress 
MPa  (pal) 


32.3  (+4690) 
31.6  (+4580) 
31.0  (+4500) 

31.4  (+4550) 


-394.0  (-57107 


+15.9  (+2300) 
+15.2  (+-2200) 
+13.1  (+1898) 
+11.2  (+1621) 
+12.6  (+1829) 


-354.2  (-51375) 
-351.9  (-50140) 
-343.0  (-49760) 
-342.4  (-49660) 
-345.1  (-50060) 


Observed  Error 
10  Tests 
MPa  (psl) 

Statistical  Error 
10  Tests 
Eq.  3.42 
MPa  (psi) 

±1.40  (±203) 
± .82  (±119) 
±1.21  (±173) 
±1.54  (±223) 

±0.97  (±141) 
±1.18  (±171) 
±1.22  (±177) 
±1.34  (±195) 

±10.04  (±1457) 
±10.98  (±1574) 
±6.89  (±1000) 
± 6.38  (±  926) 
±6.12  (±  888) 

±3.45  (±500) 
±3.96  (±574) 
±4.02  (±583) 
±4.31  (±625) 
±4.33  (±628) 

±13.70  (±1987) 
±11.44  (±1660) 
± 9.29  (±1347) 
± 6.28  (±  911) 
± 6.57  (±  953) 

H.57  (±663) 
±6.20  (±899) 
±6.71  (±973) 
±7.35  (±1066) 
±7.43  (±1077) 

±28.68  (±4160) 
±14.41  (±2090) 
±13.27  (±1925) 
±11.20  (±1624) 
±12.80  (±1857) 

±9.34  (±1325) 
±11.79  (±1710) 
±13.06  (±1894) 
±13.51  (±1959) 
±14.13  (±2049) 

±33.77  (±4898) 
±21.37  (±3099) 
±15.69  (±2276) 
±12.80  (±1856) 
±12.34  (±1791) 

±9.21  (±1336) 
±11.27  (±1635) 
±11.90  (±1726) 
±12.40  (±1799) 
±12.53  (±1818) 

±38.00  (±5513) 
±36.97  (±5362) 
±23.10  (±3351) 
±17.56  (±2547) 
±17.48  (±2536) 

±12.20  (±1770) 
±14.82  (±2150) 
±15.68  (±2275) 
±16.41  (±2380) 
±16.30  (±2365) 

* 

W rcpraaanta  the  breadth  at  half  the  maximum  intensity  in  °20. 
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counting  error.  This  phenomena  becomes  especially  pronounced  for  broad 
profiles;  e.g.,  compare  the  1045-2,  TBA  G-5  and  1090-2  samples  In 
Table  3.7. 

The  reason  for  this  result  becaasi  apparent  when  one  examines  the 
entire  procedure  for  locating  the  peak  position.  Say  the  positions 
of  the  observation  points  were  found  only  once  and  for  each  repeated 
measurement  the  data  was  taken  for  the  specified  time  at  these  same 
observation  points.  One  would  expect  the  observed  error  for  the  three 
point  fit  to  be  identical  with  the  counting  error  for  the  10  measurements. 
Indeed,  this  was  shown  to  be  the  case  in  Sec.  3.4.3. 

If,  however,  during  each  measurement  the  position  of  the  data 
points  are  rechosen,  as  is  the  case  here  because  of  the  automation 
employed  for  collecting  the  data  in  Table  3.8,  the  position  of  the  data 
points  change  due  to  random  fluctuations  in  the  observed  intensities. 

Only  if  the  actual  data  were  a perfect  parabola  would  one  expect  the 
observed  and  counting  errors  to  be  identical.  This  might  lead  one  to 
think  the  profile  swat  not  be  a perfect  parabola  in  the  region  used 
for  the  curve  fit.  It  was  seen  in  Sec.  3.3.3  that  not  subtracting  the 
background  can  lead  to  errors  in  tha  determination  of  the  region  of 
fit  of  the  parabola.  As  the  breadth  increases,  the  ratio  of  peak 
intensity  to  background  intensity  decreases  for  these  samples  (see 
Table  2.2)  causing  the  region  of  fit  as  determined  by  the  'top  15  pet' 
rule  to  be  too  large.  The  data  does  not,  then,  necessarily  lie  on  a 
parabola. 

When  a least-squares  multiple  point  fit  la  used,  the  relative 
effect  of  those  points  lying  outside  the  true  region  of  the  parabola 


I 

I 
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is  decreased  causing  Che  observed  error  Co  decrease,  alchough  Che  same 
cocal  time  of  daca  accumulation  is  used.  Surprisingly  Che  scaciscical 
councing  error  increases  probably  because  Che  degree  of  perfecdon  of 
Che  parabolic  fic  decreases  vich  increasing  number  of  daca  poinCs. 

This  is  due  Co  random  fluccuadons  in  Che  daca  which  dictate  chac  Che 
observed  daca  caimoc  exactly  lie  on  a parabolic  curve. 

To  determine  if  Che  region  of  curve  fitting  actually  effects  Che 
reproducibility,  replicate  measurements  were  made  on  che  1045>1,  1045-2 
and  TBA  G-5  samples  using  Che  same  tocal  time  as  in  Table  3.8  but  including 
a background  correction  (see  Sec.  4,3.4  for  details  on  automated  back- 
ground subtraction)  prior  Co  determining  the  region  of  fit.  The  two 
tilt  method  using  a three  point  parabolic  curve  to  define  the  peak 
location  was  used.  The  results  given  in  Table  3.9  indicate  an  improvesient 
by  a factor  of  alamst  two  in  the  observed  error  as  compared  to  Table  3.8 
for  the  samples  exhibiting  broad  profiles  when  using  a three  point 
parabolic  fit.  However,  the  precision  for  the  1045-2  and  TBA  G-5 
samplas  was  still  almost  twice  that  predicted  by  counting  statistics  alone, 
again  indicating  that  the  three  point  fit  is  not  adequate  for  fitting 
broad  profiles. 

3.6.4  Comparison  of  'Sin^f'  Method  to  'Two  Tilt'  Method 
2 

The  sin  f method  Involves  measurement  of  the  peak  position  at 
more  than  two  inclinations  which  could  mlnimsle  the  random  errors 
associated  with  individual  measurements  of  the  peak.  For  a fixed  total 
time,  the  time  spent  at  each  f tilt  will  be  less  than  in  the  two-tilt 
method,  thereby  increasing  the  statistical  error  of  the  peak  location  so 
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TABLE  3.9 

RESULTS  OF  REPLICATE  FIXED  TIME  TESTS  USING  BACKGROUND  SUBTRACTION 
(two-tilt,  three  point  perebolic  fit  to  peek) 


Staple 

Total  Tiae 
(eec) 

Street 

MPa  (pel) 

Obeerved  Error 
10  Teat 
MPa  (pai' 

Statistical  Error 
Eq.  3.42 
MPa  (pai) 

1045-1 

360 

-169.0  (-24520) 

+ 4.58  (+  605) 

+3.86  (+  560) 

1045-2 

1000 

-408.5  (-59260) 

+14.44  (+2095) 

±9.41  (+1365) 

TBA  G-5 

1000 

+11.17  (+1620) 

+17.48  (+2535) 

*9.79  (±1420) 
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that  one  cannot  say  a priori  whether  or  not  the  sin  f method  Is  more 
* 

precise. 

2 

The  precision  of  the  sin  f method  was  determined  in  the  same  manner 
as  described  for  the  two  tilt  method,  the  total  time  of  data  accumulation 
being  identical.  Only  a three  point  parabolic  fit  (no  background  sub- 
traction) was  employed,  the  ten  replicate  measurements  being  made  on 
the  1045-1,  1045-2  and  TBA  G-5  samples  with  * tilts  of  0°,  26.57°,  39.23°  and  50.77°. 

If  the  results  given  in  Table  3. 10  are  compared  to  those  in  Table 

2 

3.8,  it  is  seen  that  the  sin  f technique  has  better  precision  (observed 
error)  over  the  ten  measurements  than  the  two  tilt  technique  even  when 
the  total  time  of  data  accusailatlon  is  identical.  The  effect  of  errors 
in  peak  location  are  minimized  with  ssiltiple  tilts.  Increasing  the 
number  of  f inclinations  to  more  than  four  was  tested  but  improvement 
in  the  precision  was  only  nominal.  The  standard  error  in  the  mean 
stress  (column  4,  Table  3.10)  is  still  30  pet  greater  than  that  pre- 
dicted by  counting  statistics  (column  5)  for  the  three  point  fit  on  the 
broader  profiled  samples. 

The  best  results  when  using  a three  point  parabolic  fit  could  be 

2 

expected  to  be  found  by  combining  both  the  'sin  method  and  using  a 
background  correction  in  determining  the  region  of  peak  fit.  Ten 
replicate  measurements  were  taken  in  this  manner  on  the  same  samples 
and  using  the  sasw  fixed  time  as  in  Table  3.10.  The  results  from  such 
measurements  are  given  in  Table  3.11.  The  precision  is  indeed  improved 
and  coispares  quite  well  with  that  expected  from  counting  statistics  alone. 

* 2 
This  entire  discussion  assusms  that  the  d^  va.  sin  f data  falls  on  a 

straight  line  so  that  the  two  tilt  method  is  valid  and  accurate  . If 

df  vs.  sln*f  ware  not  linear,  the  Marl on- Cohen  method  would  have  to  be 

used  to  account  for  the  oscillations  due  to  deformation  texturing. 
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TABLE  3.10 

PRECISION  USING  SIN2*  METHOD 
(three  point  parabolic  fit  to  peak) 


Sample 

Total  Time 
(aec) 

Streaa 

MPa  (pal) 

Obaerved  Error 
10  Teata 
MPa  (pal) 

Statlatical  Error 
Eq.  3.39 
MPa  (p8i) 

1045-1 

360 

-193.0  (-27990) 

± 5.65  (±  820) 

± 5.97  (±  866) 

1045-2 

1000 

-402.8  (-58430) 

±15.93  (±2310) 

±10.90  (±1581) 

TBA  G-5 

1000 

+22.8  (+  3300) 

±18.13  (±2630) 

±14.48  (±2100) 

67 


0 


TABLE  3.11 

PRECISION  USING  SIN2*  METHOD  AND  BACKGROUND  CORRECTION 
(three  point  parabolic  fit  to  peak) 


Sample 

Total  Time 
(sec) 

Stress 
MPa  (psl) 

Observed  Error 
10  Tests 
MPa  (psl) 

Statistical  Error 
Eq.  3.39 
MPa  (psl) 

1045-1 

360 

-199.60  (-28952) 

±6.19  (±  898) 

± 4.78  (±  694) 

1045-2 

1000 

-424.20  (-61529) 

±11.17  (±1620) 

±10.61  (±1539) 

TBA  G-5 

1000 

+ 14.71  (+  2134) 

±16.13  (±2340) 

±16.65  (±2415) 

68 


i 


! 


I 


3.6.5  Precision  Obtained  Using  a Fosition  Sensitive  Detector 

In  Sec.  3.6.3  it  was  shown  that  many  data  points  within  the  region 
of  curve  fitting  improved  the  observed  precision  of  the  stress  measurement 
for  a parabolic  curve.  The  position  sensitive  detector  (PSD)  collects 
data  simultaneously  across  the  diffraction  profile  and  therefore  does 
not  take  any  longer  period  of  time  in  collecting  many  points  compared 
to  just  three.  The  number  of  data  points  collected  in  a fixed  period 
of  time  depends  solely  on  the  breadth  of  the  diffraction  peak  and  a 
least-squares  parabolic  fit  may  be  made  to  all  those  points  falling 
within  the  appropriate  region. 

Replicate  measurements  were  made  with  the  PSD  to  test  the  actual 
precision.  The  automated  residual  stress  program  was  used  without 
any  modifications.  The  fixed  time  of  analysis  was  pre-determined  by 
the  time  selected  on  the  multichannel  analyzer  (MCA)  which  acted  as  an 
Intermediate  storage  device  for  the  data  from  the  PSD.  The  MCA  was 
interfaced  with  the  PDP-8E  computer  for  direct  data  transfer. 

The  computer  program  determined  the  region  of  curve  fitting  by 
the  'top  15  pet'  of  the  intensity  profile  without  subtracting  back- 
ground and  fit  a least-squares  parabola  to  all  data  within  the  region. 
Column  2,  Table  3.12  gives  the  number  of  data  points  within  the  region 
of  curve  fitting.  The  two  tilt  method  was  utilized,  the  total  time  of 
both  tilts  being  recorded  in  column  3 of  Table  3.12. 

The  results  tabulated  in  Table  3.12  indicate  that  the  accuracy 

of  the  PSD  system  is  excellent.  The  observed  error  (column  5)  in  the 

mean  stress  is  very  close  to  that  predicted  by  the  counting  statistics 

2 

(column  6)  alone.  Thus  the  sin  f method  is  not  needed  if  the  'two  tilt' 
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TABLE  3. 12 


PRECISION  OBTAINED  WITH  A POSITION  SENSITIVE  DETECTOR 
(TVo-Tilt  Method) 


Sample 

Number  of 
Data  Polnta 

Time 

(aec) 

S trees 
MPa  (psi) 

Observed  Error 
10  Tests 
MPa  (psi) 

Statistical  Error 
MPa  (psi) 

1090-1 

15 

100 

23.89  (+  3465) 

6.34  (±  920) 

4.48  (±  650) 

1045-1 

31 

100 

198.6  (-28800) 

5.72  (±  830) 

4.34  (±  630) 

1045-3 

87 

200 

718.4  (-104200) 

7.38  (±1070) 

6.87  (±  995) 

1045-2 

141 

200 

412.3  (-59800) 

9.20  (±1335) 

9.69  (±1405) 

TBA  G-5 

196 

200 

22.75  (+  3300) 

13.03  (±1890) 

13.86  (±2010) 
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method  is  accurate  (see  Sec.  1.2).  Obtaining  many  data  points  (column  2) 
within  the  region  of  parabolic  fit  eliminates  the  necessity  for  sub- 
tracting background.  Hie  best  precision  obtained  with  the  normal 

2 

detector  comes  from  using  the  sin  f method  and  subtracting  background. 

The  time  for  the  final  data  collection  for  this  procedure  (not  Including 
the  60  seconds  or  so  taken  during  the  preliminary  scans)  is  given  in 
Table  3.11  and  can  be  compared  with  the  total  time  of  analysis  using 
the  PSD  in  Table  3.12.  This  demonstrates  the  remarkable  speed  of  the 
PSD,  even  when  only  a normal  X-ray  tube  is  used. 

3.7  STUDIES  OF  SAMPLE  DISPLACEMENT 

Perhaps  the  largest  single  source  of  instrumental  error  aasociated 

with  the  X-ray  measuresmnt  of  residual  stresses  is  sample  displacesttnt. 

In  focusing  geometry,  if  the  effective  diffracting  volume  is  not  located 

at  the  center  of  the  diffraction  circle,  there  is  a relative  shift 

between  f-0°  and  given  by  Eq.  3.6.  The  error  due  to  a given 

displacement  la  dependent  on  R the  poaltlon  of  the  receiving  slit. 

P 

The  equation  is  valid  for  both  the  parafocualng  and  stationary-slit 
techniques  and  is  experimentally  tested  below.  The  parallel  beam 
technique  discussed  in  Sec.  3.2.2  does  not  employ  a receiving  slit  to 
define  the  angular  relationship  between  the  primary  and  diffracted 
beam.  The  angular  relationship  depends  solely  on  the  angle  between  the 
parallel  Soller  slits  nwunted  90°  to  the  diffractometer  plane  and  does 
not  depend  on  the  position  of  the  sample,  thereby  eliminating  the  error 
caused  by  sample  dlsplacesient. 

The  effect  of  saaple  displacement  was  tested  for  each  of  the 
"beam  optics"  techniques.  The  positioning  of  the  sample  was  accurate 
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Co  within  ± .0025  ran.  Three  replicate  measurements  were  made  at  each 
sample  position  up  to  a total  displacement  of  i 2 ran*  in  increments  of 
.5  ran  using  the  automated  residual  stress  program.  The  1045-1  sample 
was  utilized  because  it  gave  a reasonably  sharp  profile.  The  residual 
stress  was  measured  using  the  two  tilt  procedure  and  a three  point 
parabolic  fit  without  background  correction. 

The  parallel  beam  method  was  obtained  by  placing  two  sets  of  high  reslution 
Picker  Soller  slits  together  in  both  the  incident  aad  diffracted  beams  and 
rotated  90°  to  the  usual  position.  This  yielded  a divergence  of  .5° 
similar  to  the  values  used  in  Japan.  Such  slit  systems  were  placed 

in  the  primary  and  diffracted  beams. 

The  error  for  both  the  parafocusing  and  stationary  slit  methods 
predicted  by  Eq.  3.6  are  drawn  as  solid  lines  in  Fig.  3.7.  The  experi- 
mental data  given  In  Fig.  3.7  for  each  technique  closely  follows  the 
predicted  error.  The  stationary  slit  technique  is  seen  to  be  less 
susceptible  to  sample  displacement  by  almost  a factor  of  five  over  the 
parafocusing  technique. 

The  parallel  beam  geometry  la  seen  to  be  insensitive  to  reasonable 
sample  displacements  in  agreement  with  Aoyama , et  al.*68)  fof  course, 
if  the  sample  is  displaced  far  enough,  the  diffracted  beam  will  not  fall 
completely  in  the  receiving  Soller  slits  introducing  a large  error.] 

The  parallel  beam  technique  decreases  the  intensities  by  about  % and 
broadens  the  diffraction  profile.  This  instrumental  broadening  is  oust 
important  on  samples  exhibiting  sharp  profiles  when  focusing  geometry  is  used. 
The  reproducibility  or  precision  of  the  parallel  bean  technique  was  tested. 

* 

A negative  displacesMnt  is  defined  as  being  towards  the  back  surface  of 

the  sample. 


The  data  for  Table  3.13  was  acquired  for  a fixed  time  to  give  a statis- 
tical error  close  to  that  in  Table  3.8  using  the  two  tilt  procedure  and 
a three  point  fit  (no  subtraction  of  background).  The  precision  for 
the  1045-1  is  twice  as  poor  but  approaches  that  of  the  focusing  techniques 
for  samples  having  broad  profiles.  The  time  of  analysis  is  1/3  longer  for 
approximately  the  same  statistical  error  for  all  specimens. 

3.8  SUMMARY 

A review  of  the  important  biasing  factors  in  the  residual  stress 
measurement  by  X-ray  diffraction  techniques  has  been  presented.  The 
effect  of  sample  displacement  has  been  treated  experimentally  verifying 
the  formulae  used  to  predict  this  error.  It  was  seen  that  stationary 
slit  'beam  optics'  is  less  susceptible  to  sample  displacement  than 
parafocuslng  geometry  and  has  leas  Instrumental  broadening  than  the 
parallel  beam  method.  This  indicates  the  stationary  slit  method  la  the 
best  'beam  optics'  arrangement  when  applied  to  X-ray  residual  stress 
aeasureswnts . 

A complete  statistical  analysis  for  the  error  in  stress  due  to 

counting  statistics  was  presented.  The  analysis  was  made  as  general  as 

possible  and  includes  unlimited  data  points  in  the  peak  region  and  use 

2 

of  either  two  tilt  or  'sin  f'  procedures. 

It  has  been  shown  experimentally  that  the  three  point  parabola 
(which  la  traditionally  used  because  of  its  simplicity)  la  justified 
only  when  the  region  of  curve  fit  la  determined  accurately.  If  the 
background  level  la  high  compared  with  the  peak  intensity.  It  must  be 
subtracted  prior  to  determining  the  region  of  fit  using  the  'top  15  pet' 
rule.  Even  If  this  la  accomplished,  a three  point  fit  la  not  very 
precise  when  broad  proftles  are  Involved.  By  using  multiple  data  points 
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TABLE  3.13 


PRECISION  USING  PARALLEL  BEAM  GEOMETRY 
(TWo  Tilt  Procedure  - Three  Point  Parabolic  Fit  to  Peak) 

(10  Measuremants) 


Sample 

Bread  th 
(°20) 

Time 

(aec) 

S trees 

MPa  (pal) 

1045-1 

2.1 

500 

-161.0  (-  23360) i 

1045-3 

3.8 

1300 

-702.6  (-101920) j 

1045-2 

5.3 

1300 

-417.3  (-  60525) 

TBA  G-5 

HI 

1300 

-1-  11.6  (+  1682) 

Statistical  Error 
MPa  (pal) 
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and  • least-squares  parabola  on  broad  profiles  the  reproducibility  of 
peak  location  can  be  improved  by  a factor  of  three  without  increasing 
the  total  time  of  data  accumulation.  The  practical  necessity  of  using 
a computer  to  perform  the  least-squares  analysis  detracts  from  using 
multiple  data  points,  however,  when  complete  autoaution  is  available 
this  type  of  analysis  is  worthwhile.  In  addition,  the  practical  advan- 
tage of  multiple  data  points  increases  the  appeal  of  a position  sensitive 
detector  where  such  data  is  accumulated  simultaneously. 

2 

Surprisingly,  when  only  a three  point  fit  is  used,  the  sin  f 

method  is  shown  to  be  more  reproducible  than  the  two  tilt  method  even 

if  the  total  time  is  fixed  in  both  procedures.  Clearly,  if  true 

automation  is  to  be  achieved  with  optimum  procedures,  a least-squares 

2 

parabola  to  define  the  peak  or  the  sin  f method  of  analysis  must  be 
used. 


D 


76 


REFERENCES 


1.  Nondestructive  Evaluation  of  Residual  Stress,  proceedings 
of  s Workshop,  Air  Force  Materials  Laboratory,  NTIAC-76-2, 
Southwest  Research  Center,  San  Antonio,  Texas  (1976). 

2.  B.  V.  Weinman,  J.  B.  Hunter,  and  D.  D.  McCormack,  Metal  Progress, 
96  No.  1,  88  (1969). 

3.  K.  Ramachl  , X-Ray  Study  on  Strength  and  Deforsation  of  Metals, 
The  Society  of  Materials  Science,  Japan  (1971),  p.  95. 

4.  C.  S.  Barrett  and  T.  B.  Massalskl,  Structure  of  Metalg,  3rd 
Edition,  McGraw  Hill,  New  York  (1966),  p.  466-485. 

5.  B.  D.  Cullity,  Elements  of  X-Ray  Diffraction,  Addison-Uesley 
Publ . Co.,  Reading,  MA.  (1956),  p.  431-451. 

6.  H.  P.  Klug  and  L.  E.  Alexander,  X-Ray  Diffraction  Procedures, 
John  Wiley  and  Sons,  Inc.,  New  York  (1974)  p.  755-790. 

7.  Residual  Stress  Measurement  by  X-Ray  Diffraction-SAE  J 784a, 

2nd  Edition,  Society  of  Automotive  Engineers,  Inc.  (1971). 

8.  6.  B.  Greenough,  Prog,  in  Metal  Physics,  3,  176  (1952). 

9.  S.  Taira,  R.  Hayashi  and  Z.  Watase,  Proceedings  of  the  IVelfth 
Japan  Congress  on  Materials  Research,  The  Society  of  Materials 
Science,  Japan,  1 (1970). 

10.  T.  Shiralwa  and  Y.  Sakamoto,  Proceedings  of  the  Thirteenth 
Japan  Congress  on  Materials  Research,  The  Society  of  Materials 
Science,  Japan,  25  (1971). 

11.  E.  Macherauch,  Proc.  3rd.  Int.  Conf.  Nondes.  Test.,  720  (1961) 
Tokyo,  Japan. 

12.  A.  Ruess,  Z.  Angew,  Math.  Mech.,  9,  p.  48-58  (1929). 

13.  W.  Volght,  Lehrbuch  der  Krlstallphyslk,  Teubner,  Leipzig/Berlin, 
(1928) . 

14.  B.  Kroner,  Z.  Phys.,  151,  p.  504-508,  (1958). 


r 


77 


15.  F.  Bollenrath,  V.  Hauk,  and  E.  H.  Muller,  Z.  Metallkunde,  58, 

p.  76-82  (1967).  ~ 

16.  N.  Hosokawa,  K.  Honda,  and  J.  Arlma,  Mechanical  Behavior  of 
Materials,  (Proceedings  of  the  Inter.  Conference),  The  Society 
of  Materials  Science,  Japan,  p.  164-175  (1972),  Kyoto,  Japan. 

17.  J.  T.  Morton,  Adv.  In  X-Ray  Analysis.  11,  p.  401  (1967). 

18.  R.  E.  Ricklefs  and  W.  P.  Evans,  Adv.  in  X-Ray  Analysis, 

10,  p.  273  (1966). 

19.  R.  H.  Marion  and  J.  B.  Cohen,  Adv.  in  X-Ray  Analysis,  IB, 
p.  466  (1974). 

20.  T.  Shlraiva  and  Y.  Sakamoto,  X-Ray  Study  on  Strength  and  De- 
formation of  Metals,  The  Society  of  Materials  Science,  Japan, 

(1971)  p.  15. 

21.  F.  Bollenrath,  V.  Hauk,  and  W.  Ohly,  Z.  Metallkunde,  57, 
p.  464  (1966). 

22.  M.  J.  Donachie  and  J.  T.  Horton,  Trans  AIMS,  221,  p.  962 
(1961). 

23.  F.  Bollenrath,  V.  Hauk  and  U.  Weldemann,  Archiv  fur  das 
Eisenhtftt,  38  793  (1967). 

24.  V.  Weldemann,  Fh.D.  thesis,  Technische  Hochschule,  Aachen, 

Germany  (1966) . 

25.  K.  W.  Andrews,  J.  C.  Gregory  and  D.  Brooksbank,  Strain,  10, 

111  (1974). 

26.  D.  P.  Koistinen  and  R.  E.  Marburger,  Trans  ASM,  51,  537  (1959). 

27.  X-Ray  Studies  on  Mechanical  Behavior  of  Metals,  The  Society 
of  Materials  Science,  Japan,  Kyoto,  Japan  (1974)  p.  72. 

28.  M.  Short  and  C.  Kelly,  Adv.  In  X-Ray  Analysis,  16.  379  (1972). 

29.  C.  Jatczak  and  H.  Boehm,  Adv.  in  X-Ray  Analysis.  17,  354  (1973). 

30.  R.  H.  Marlon,  Fh.D.  thesis.  Northwestern  University,  1973. 

31.  C.  Jatczak  and  H.  Zantopulos,  Adv.  1>,  X-Ray  Analysis.  14,  360  (1970) 

32.  Residual  Stress  Measurements  by  X-Ray  Diffraction  - SAE  J784a, 

2nd  Edition,  Society  of  Automotive  Engineers,  Inc.  (1971)  p.  71. 

J»  B.  Cohen,  Diffraction  Methods  in  Materials  Science,  MacMillan  Co. 
Mew  York  (1966)  p.  294. 


33. 


78 


34.  K.  H.  Lauterjung,  J.  Pokar,  B.  Schlmmer,  and  R.  Staudner, 

Nucl.  Inst,  and  Meth.,  22,  117  (1963). 

35.  R.  Bock,  H.  H.  Duhn,  U.  Melzer,  F.  Dulilhofer,  and  B.  Stadler, 

Nucl.  Inst,  and  Meth.,  41,  190  (1966). 

36.  S.  Kalbltzer  and  W.  Melzer,  Nucl.  Inst,  and  Meth.,  56,30  (1967). 

37.  W.  W.  Daehnick,  Phys.  Rev.,  177,  1795  (1967). 

38.  D.  L.  Alan,  G.  V.  Arsell,  and  R.  K.  Kelly,  Atomic  Energy  Research 
Est.,  Harwell,  Berkshire  Report  AERE-R5897  (1968). 

39.  A.  Doehrlng,  S.  Kalbltzer,  and  W.  Melzer,  Nucl.  Inst,  and  Meth. 

59,  40  (1968) . 

40.  R.  B.  Owen  and  M.  L.  Awcock,  IEEE  Trans.  Nucl.  Scl.,  NS-15,  290 
(1968). 

41.  U.  Beble  and  F.  Parah,  Nucl.  Inst,  and  Meth.,  112,  455  (1973). 

42.  W.  R.  Kuhlmann,  K.  H.  Lauterjung,  B.  Schlmmer,  and  L.  Slstermlch, 

Nucl.  Inst,  and  Meth.,  40,  118  (1966). 

43.  U.  N.  McDlchen,  Nucl.  Inst,  and  Meth.,  54,  157  (1967). 

44.  C.  J.  Borkowskl  and  M.  K.  Kopp,  Rev.  Scl.  Inst.,  39,  1515  (1968). 

45.  C.  J.  Borkowskl  and  M.  K.  Kopp,  IEEE  Trans.  Nucl.  Scl.,  L7,  340  (1970). 

46.  J.  Hough  and  A.  Drever,  Nucl.  Inst,  and  Meth.,  103,  365  (1973). 

47.  H.  Foeth,  R.  Hammer strom,  and  C.  Rubbla,  Nucl,  Inst,  and  Meth., 

109.  521  (1973). 

48.  D.  G.  Smith  and  K.  A.  Pounds,  IEEE  Trans.  NUcl.  Scl.,  15,  541  (1968). 

49.  R.  Gott,  W.  Parks  and  K.  A.  Pounds,  Nucl.  Inst,  and  Meth.,  81,  152 

(1970) . “ 

50.  J.  Hough,  Nucl.  Inst,  and  Meth.,  105.  323  (1972). 

51.  G.  I.  Miller,  A.  Senator,  and  R.  Stensgraad,  Nucl.  Inst,  and  Meth., 

91,  389  (1971). 

52.  G.  P.  Westphal,  Nucl.  Inst,  and  Meth.,  106.  279  (1973). 

53.  H.  W.  Pulbrlght,  R.  G.  Markham,  and  W.  A.  Langford,  Nucl.  Inst,  and 
Meth.,  108,  125  (1973). 

54.  J.  L.  C.  Ford,  J.  Gomes,  Del  Caaqpo,  R.  L.  Robinson,  B.  H.  Stelson, 
and  S.  T.  Thornton,  Nuclear  Physics,  A226,  189  (1974). 


55.  A.  Gabriel  and  S.  Bram,  Febs.  Letters,  39,  307  (1974) 


Y.  Dupont,  A.  Gabriel,  MChabre,  T.  Gulik-Krzyvlckl,  and 
E.  Schechter,  Nature,  238,  331  (1972). 


57.  E.  Mathleson,  Nucl.  Inst,  and  Meth.,  97,  171  (1971) 


H.  P.  Klug  and  L.  E.  Alexander,  X-Ray  Diffraction  Procedures 
John  Wiley  and  Sons,  Inc.,  New  York  (1974)  p.  329. 


60.  F.  Y.  Yap,  dissertation  (The  John  Hopkins  University 
Baltimore,  Maryland,  1967). 


61.  A.  J.  C.  Wilson,  Acta.  Cryst.  22,  888  (1967) 


J.  Thomsen  and  F.  Y.  Yap,  J.  of  Research  of  Nat.  Bureau  of 
Standards,  72A,  187  (1968). 


63.  A.  J.  C.  Wilson,  Brit.  J.  Apply.  Phys.,  16,  665  (1965). 

64.  M.  J.  Cooper  and  A.  V.  Glasspool,  J.  Apply.  Cryst.  9,  63  (1976) 


International  Tables  for  X-Ray  Crystal lograph 
Kynoch  Press,  Birmingham,  England  (1968). 


J.  B.  Cohen,  "Report  on  Tungsten  Lattice  Parameter  Round  Robin," 
X-Ray  Subcommittee  of  SAE  Iron  and  Steel  Technical  Committee, 
Division  4,  (November  1964). 


J.  Fukura  and  H.  Fujlwara,  J.  of  the  Society  of  Materials  Science 
Japan,  15,  825  (1966). 


S.  Aoyama,  K.  Satta  and  M.  Tada,  J.  of  the  Society  of  Materials 
Science,  Japan,  17,  1071  (1968). 


69.  D.  N.  French,  J.  Amer.  Cer.  Soc.,  52,  271  (1969) 


W.  E.  Baucum  and  A.  M.  Ammons,  Adv . in  X-Ray  Analysis  .17.371 
(1973'.. 

"Standard  Method  for  X-Ray  Stress  Measurement",  The  Committee  on 
Mechanical  Behavior  of  Materials,  The  Society  of  Materials  Science 
Japan  (1973). 


72.  Residual  Stress  Maasuri 
2nd  Edition,  Society  o: 


>nts  by  X-Ray  Dlffractlon-SAE  J784a, 
automotive  Engineers,  Inc.,  (1971),  p.  51 


Science,  Japan  (l1 


80 


74.  A.  J.  C.  Wilson,  Mathematical  Theory  of  X-Ray  Powder  Dlffracto- 
metry , Philips  Technical  Library  (1963) . 

75.  W.  Parrish  and  T.  C.  Huang,  Profile  Tilting:  A Powerful  Method 

of  Computor  X-Ray  Instrumentation  and  Analysis,  IBM  Research  (1976). 

76.  0.  Davies  and  P.  Goldsmith,  Statistical  Methods  In  Research 
and  Production,  Hefner  Publ.  Co.,  Hew  York  (1952)  p.  72. 

77.  X-Ray  Diffraction,  Ed.  L.  Araroff,  McGraw  Hill,  New  York  (1974) 
p.  439. 


DOCUMENT  CONTROL  DATA  ■ R & D 


(Sfrurify  clmaatllcation  ol  title,  body  of  «b»fr*r^ii^nrfeji/H£Annof*»|o^ 

" |2«.  REPORT  SECUrTT^ 


Unclassified 


J.  B.  Cohen,  Northwestern  University,  Evanston,  IL 


EVALUATION  OF  PROCEDURES  IN  AUTOMATED  RESIDUAL  STRESS  MEASUREMENTS 


« ocicriptive  noth  (7>p#  of  report  o nO  Inctumtv * drnij 

Technical  Report  No.  14 

• kUTHOmil  enrol  nMM,  aitaldfo  Initial,  foot  namo) 


M.  R.  James  and  J.  B.  Cohen 


N00014-75-C-0580 


Technical  Report  No.  14 


OTHER  REPORT  NOIJI  (Any  other  number  a that  may  be  meelgned 
thte  report) 


None 


Distribution  of  this  document  is  unlimited 


Office  of  Naval  Research  Metallurgy  Branch 


^^^'"The  measurement  of  residual  stress  by  the  X-ray  diffraction  technique  involves 
determining  the  lattice  spacing  of  a crystallographic  plane  at  different  inclinations 
of  the  sample  and  relating  the  change  in  the  spacing  to  a stress  on  the  surface  of 
the  specimen.  The  sources  of  error  in  determining  the  residual  stress  are  inves- 
tigated in  this  report.  A short  review  of  the  fundamental  principles  of  the 
measurement  is  given.  The  important  Instrumental  and  geometric  factors  contributing 
to  errors  in  the  measured  stress  are  presented.  To  account  for  random  errors  in  the 
data  accumulation  on  the  measured  stress,  a complete  statistical  analysis  based  on  a 
least-squares  parabola  is  given. 

tome ter  (see  T.R 


Employing  an  automated  dlff  fl^tometer  (see  T.R.  No.  16),  an  extensive  investi- 
gation of  the  precision  of  the  r^slan^l  stress  measurement  was  made  to  evaluate  the 
various  procedures  offered  in  the  literature  to  optimize  the  automation.  The  results 
indicate  that  a multiple  least-squares  parabola  is  the  most  reproducible  method  of 
defining  the  peak  of  the  profile,  in  contrast  to  the  common  three  point  parabola 
used  in  the  United  States,  the  center  of  gravity  used  in  Germany  or  the  middle  of  the 
half-maximum  intensity  used  in  Japan.  Also,  for  a standard  diffractometer,  the  s In® 
technique  la  more  precise  than  the  usual  two  tilt  method  even  when  the  total  time  of 
analysis  for  each  is  identical.  In  studying  various  beam  optic  arrangements,  it  was 
found  that  the  stationary  slit  (non-focusing  method)  offers  the  best  system  in  that 
the  omasured  stress  is  only  mildly  sensitive  to  sample  displacement  (~  7 MPa  for  .25 
mm  displacement)  and  introduces  very  little  Instrumental  broadening.  The  operational 
principles  and  use  of  a position  sensitive  detector  (PSD)  for  the  measurement  of  resi 
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